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We discuss first some basic experimental facts related to ARPES, tunnelling, optics ad neutron 
scattering measurements. They give evidence for the relevance of the electron-phonon interaction 
(EPI) in pairing mechanism of HTSC cuprates. A controllable and very efficient theory for strong 
correlations and their effects on EPI is discussed. The theory is based on the 1/N expansion 
method in the t-J model without using slave bosons (or fermions). The remarkable prediction 
is that strong correlations renormalize EPI and other charge-fluctuation properties (by including 
nonmagnetic impurity scattering) in such a way that the forward scattering peak (FSP) appears 
in the corresponding effective interactions. The existence of FSP in EPI is additionally supported 
by the weakly screened Madelung interaction in the ionic-metallic structure of layered cuprates. 
Pronounced FSP in EPI of HTSC cuprates reconciles many puzzling results, which could not be 
explained by the old theory with the momentum independent EPI. For instance, EPI with FSP gives 
that the couplings in the s- and d-wave pairing channel are of the same magnitude near and below 
the optimal hole doping. It is shown that FSP in the impurity scattering potential is responsible 
for robustness of d-wave pairing in cuprates with nonmagnetic impurities and other defects. The 
ARPES kink and the isotope effect in the nodal and anti-nodal points, as well as the collapse of the 
elastic nonmagnetic impurity scattering in the superconducting state, are explained by this theory 
in a consistent way. The proposed theory also explains why the nodal kink is not-shifted in the 
superconducting state while the anti-nodal kink is shifted by the maximal superconducting gap. It 
turns out that in systems with FSP in EPI besides the classical phase fluctuations there are also 
internal fluctuations of Cooper pairs. The latter effect is pronounced in systems with long-ranged 
pairing forces, thus giving rise to an additional contribution to the pseudogap behavior. 

PACS numbers; 



I. INTRODUCTION 

In this year the physics community celebrate twenty 
years from the remarkable discovery of high-temperature 
superconductors (HTSC) by G. Bednorz and K. A. 
MiiUer. However, all of us are aware of the fact that 
our friend and teacher Vitalii Lazarevich Ginzburg P| is 
the pioneer in this fascinating field. He was not only the 
first who raised this question but he also actively advo- 
cated the possibility of the room-temperature supercon- 
ductors (RTSC) 0. Moreover, he has proposed the well 
known excitonic mechanism of pairing in quasi 2-D struc- 
tures, which, in fact, was a precursor for other studies of 
this important subject. Finally, HTSC were discovered 
in 1986 by Bednorz and MiiUer in rather exotic quasi- 
2D HTSC materials - cuprates A number of inter- 
esting approaches for reaching high critical temperature 
Tc were elaborated in the famous VL superconductivity- 
group in the I. E. Tamm department at the P. N. Lebe- 
dcv Institute in Moscow. (VL is the popular name of V. 
L. Ginzburg - see Q.) One of the most important re- 
sults, which came from the "theoretical kitchen" of the 
VL group, was related to a possible limit of Tc due to 
the electron-phonon interaction (EPI). 

Regrettably, study of the pairing mechanism in cuprate 
superconductors was significantly influenced by preju- 
dices related to EPI. One prejudice in the search of 
HTSC materials is related to the question of the up- 



er limit of Tc within the phonon mechanism of pairing 
By assuming that the total macroscopic longitudi- 
nal dielectric function StotO^,^ = 0) is positive for all k 
(£tot(k, w = 0) > , k is the momentum and uj is the 
frequency) the authors in Ref. Q came to a pessimistic 
conclusion that EPI is able to produce only a moderate 
Tc < 10 K, since as they predicted the EPI coupling con- 
stant is limited to its maximal value X']^^j < 0.5. How- 
ever, a number of experiments contradicted this state- 
ment, since the large coupling constant Xepi > 1 is real- 
ized in a number of materials. For instance Xepi ~ 2.5 in 
the PbBi alloy (although Tc is not large) , and A epi ~ 1 
with Tc 4:0 K in MgB2- There are a number of other 
systems with Xepi > 0.5, such as A3 Ceo with Xepi ~ 1 
and Tc ~ 40 K, etc. The correct answer to this question 
was given by D. A. Kirzhnitz 0, a remarkable man and 
physicist from the Ginzburg group, who has shown that 
thermodynamic and lattice stability do not exclude the 
possibility that etot(k, = 0) < 0. In fact a negative 
value of etot(k, w = 0) is realized in a large portion of the 
Brillouin zone in most superconductors 6J. This impor- 
tant result means that the dielectric function Stot does 
not limit Tc in the phonon mechanism of pairing - see 
more in 0, 0, and 0. 

Concerning HTSC in cuprates one of the central ques- 
tions is - what is the pairing mechanism in these mate- 
rials? In the last twenty years, the scientific community 
was overwhelmed by all kinds of (im)possible proposed 
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pairing mechanisms, most of which are hardly verifiable, 
if at all. This trend is still continuing nowadays (although 
with smaller slope), in spite of the fact that the accu- 
mulated experimental results eliminate all but few. A 
somehow similar situation happened also in experimen- 
tal physics of HTSC cuprates, where a whole spectrum of 
possible and impossible values of some (for pairing mech- 
anism) relevant quantities was reported. As illustration 
the reader can look at (critical) analysis of experimental 
situation in optics done by Ivan Bozovic - especially of 
the measurements of the plasma frequency cOpi . The dis- 
persion in the reported value of ujpi, from 0.06 to 25 eV - 
almost three orders of magnitude - tells us that in some 
periods the science suffers from the lack of rigorousness 
and objectiveness. 

By the end of the twenty-years era in studying HTSC 
materials, more or less three main directions in search for 
the pairing mechanism have crystallized. The first ap- 
proach, based on strong correlations, was first proposed 
by P. W. Anderson in 1987 and followed by many others. 
The starting point is that the physics of cuprates is the 
physics of doped Mott insulator 10], which is supposed to 
be well described by the single-band Hubbard model (or 
its extensions). The promoters of this approach believe 
that strong electronic correlations alone are the driving 
mechanism behind the whole phase diagram of these ma- 
terials - the antiferromagnetic phase, metallic and super- 
conducting one. Since strong correlations belong to the 
strong coupling problem it is understandable that even 
nowadays, exactly 20 years after the discovery of HTSC, 
there is still no reliable theory of strong correlations in 
HTSC materials which is capable to describe both the 
normal and the superconducting state. 

The second approach is based on pronounced antiferro- 
magnetic spin fluctuations which are present in cuprates 
especially at very low doping - we call it the spin- 
fluctuation interaction approach (SFI) Later on we 
will argue that this phenomenological approach - which 
can be theoretically justified in the very weak coupling 
limit only - is based on two ingredients which are more in 
conflict than in agreement with the existing experimental 
results. First, the spin-fluctuation spectrum taken from 
the (low-frequency) NMR measurements differs from the 
one obtained in neutron scattering measurements - see 
0, 01- The latter is shifted to much higher energies. 
Second, in order to explain large Tc within the SFI ap- 
proach, one has to assume an unrealistically large spin- 
fluctuation coupling constant gsf ~ 0.6 eV, which implies 
that, the pairing coupling constant is also unrealistically 
large As/(~ g^^) > 2. Our statement is confirmed by a 
number arguments. Let us mention here some of them, 
while others are given in Sec. II. A. For instance, such a 
large Xgf gives much larger resistivity than the measured 
one. Also, gsf extracted from ARPES and neutron scat- 
tering measurements is much smaller, i.e. gsf ~ 0.2 eV 
what gives rather small pairing constant Xsf < 0.3, see 
Sec. II A. 

Let us stress that both of the above approaches are 



based either on the Hubbard model or on its (more pop- 
ular) derivative, the t-J model. Here, the central question 
is - do these two models show superconductivity at high 
temperature? The Monte-Carlo calculations on the 2D 
repulsive {U > 0) Hubbard model give no evidence for 
superconductivity with a large critical temperature |l3j |. 
On the other hand, there is a strong tendency for super- 
conductivity (either BCS like or Berezinskii-Kosterlitz- 
Touless like in 2D systems) in the attractive {U < 0) 
Hubbard model for the same intensity of U. Recent nu- 
merical calculations in the t-J model have shown that 
there is no superconductivity in the t-J model at temper- 
atures characteristic for cuprates. If it exists (at all!) 
Tc must be very low ^ I K. These numerical results 
tells us that the lack of HTSC in the repulsive {U > 0) 
single-band Hubbard model, as well as in the t-J model, 
is not due to 2D-fluctuations (which at finite T suppress 
superconductivity with phase coherence) but it is due 
to the inherent ineffectiveness of strong correlations and 
spin fluctuations to produce HTSC in cuprates. This 
means also that the simple single-band Hubbard and the 
t-J model are insufficient to explain the pairing mech- 
anism in cuprates, and some other ingredients must be 
included. Since EPI is strong in oxides it is physically 
plausible and justified that it should be accounted for. 
As we shall argue later on, the experimental support for 
the importance of EPI comes from optics, tunnelling, and 
especially from recent ARPES measurements. 

Regarding the EPI one can ask whether it contributes 
significantly to d-wave pairing in cuprates? Surprisingly, 
in spite of a number of experiments in favor of EPI many 
of researchers believe that EPI is irrelevant and that 
the pairing mechanism is due to SFI and strong cor- 
relations alone - see This belief is mainly based 
on: (i) an incorrect stability criterion, discussed briefiy 
above, and (ii) the fact that in cuprates d-wave pair- 
ing (A(k, w) « As{k,LLj) -t- Ad(w)(cosfca; — cos ky) with 
As < 0.1 Ad) is realized with gapless regions on the Fermi 
surface which is incompatible with EPI. So, the EPI 
mechanism of pairing in cuprates has been abandoned 
mainly because of these two (incorrect) prejudices. 

In fact, there is rich experimental evidence that EPI is 
sufficiently large in cuprates in order to produce HTSC, 
i.e. the EPI coupling constant is A > 1. Let us 
quote some of them: (1) The superconductivity-induced 
phonon renormalization is much larger in cuprates 
than in LTSC superconductors. This is partially due to 
the larger value of A/E^f in HTSC than in LTSC. (2) 
The line-shape in the phonon Raman scattering is very 
asymmetric (Fano line), which points to a substantial 
interaction of lattice vibrations with some quasiparticle 
continuum (electronic liquid). For instance, the phonon 
Raman measurements Q on H gBa2Ca^CunOiQj^x at 
T < Tc give very large softening (self-energy effects) of 
the Aig phonons with frequencies 240 and 390 cm~^ by 
6 % and 18 %, respectively. At the same time there 
is a dramatic increase of the line-width immediately be- 
low Tc, while above Tc the line-shape is strongly asym- 
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metric. A substantial phonon renormalization was ob- 
tained in {Cu,C)Ba2Ca^CuiOiQj^x 0- (3) The large 
isotope coefficients {ao > 0.4) in YBCO away from the 
optimal doping and ao ~ 0.15 — 0.2 in the opti- 
mally doped Lai.s^SrQ i^CuOi. At the same time one 
has ao ~ acu making a « 0.25 — 0.3. This result tell 
us that other, besides O, ionic vibrations participate in 
pairing. (4) Very important evidence that EPI plays an 
important role in pairing comes from tunnelling spectra 
in cuprates, where the phonon-related features have been 
clearly seen in the I — V characteristics [13. (5) The pen- 
etration depth in the a-b plane of YBCO is significantly 
increased by the substitution O^^ O^^, i.e. (AAafc/ 
Xab) = (l^Aah Xab)/^^\ab = 2.8 % at 4 if Sincc 
Aafe ~ fn* , the latter result, if confirmed, might be due 
to a nonadiabatic increase of the effective mass m*. (6) 
The breakthrough in the physics of cuprates came from 
the recent ARPES measurements 1131, [I^, HI. The lat- 
ter show a number of features typical for systems with 
pronounced EPI and peculiar Coulomb interaction. For 
instance, there is a kink in the quasiparticle spectrum 
in the nodal point at characteristic (oxygen) phonon fre- 
quencies in the normal and superconducting state. The 
ARPES measurements indicate that the EPI coupling 
constant is of the order one, A w Xepi ^ 1 in the nodal 
direction, and X ppj 2 near the antinodal point. The 
authors of Ref. |23| even claim that they are able to re- 
solve the fine structure in i?el](k, oj) in ARPES, what re- 
sults in the Eliashberg spectral function with a number 
of peaks at the characteristic phonon frequencies. The 
imaginary part of the ARPES self energy, ImT,{k,uj), 
has a knee-like structure around cuph ~ 70 meV. which is 
typical for systems with pronounced EPI [23|, HJ, [23l| . 
Here, LUph is the characteristic maximal phonon energy. 
For Lu > ujph one has ImYi(k,uj) ~ — AcW, which is cer- 
tainly due to Coulomb interaction (by including SFI too) 
with rather moderate Ac ^ 0.4. Recent ARPES spectra 
in the anti-nodal point give evidence for the kink at 40 
meV, while the spectral function shows a peak-dip-hump 
structure which is characteristic for systems with strong 
EPI {Xepi > 1). Very recent ARPES spectra 24] show a 
pronounced isotope effect in the real part of the quasipar- 
ticle self-energy. These results give clear evidence that 
the EPI is appreciate and strongly involved in pairing 
and quasiparticle scattering. 

On the theoretical side there are self-consistent LDA 
band-structure calculations which (in spite of their short- 
comings) give a rather large bare EPI coupling constant 
A ~ 1.5 in Lai ssSrQ i^CuOi ^2E\. The tight-binding 
calculations of EPI in YBCO gave A « 2 and « 90 
K [26j |. The nonadiabatic effects due to poor metallic 
screening along the c-axis, not accounted for in p5j |. may 
additionally increase A as pointed first in [l^. All these 
facts are in favor of a substantial EPI in cuprates. How- 
ever, if the properties of the normal and superconducting 
state in cuprates are interpreted in terms of the standard 
EPI theory - which is well established in low-temperature 
superconductors (LTSC) - some puzzles arise. One of 



them is related to the normal-state resistivity in opti- 
mally doped cuprates, which indicates that the transport 
coupling constant Xtr is rather small, i.e. At^ < A/3 ac- 
cording to ARPES and tunnelling measurements Q , ■ 
For instance, the combined resistivity and low frequency 
conductivity (Drude part) measurements give Xtr ~ 0.3 
if the plasma frequency takes the value ujpi ^ 1 eV - see 
more below. So, if one assumes that Xtr ~ A, which is the 
case in most LTSC, such a small A can not give large Tc 
(« 100 K). This implies two possibilities: (a) Xtr ^ A 
and the pairing is due to EPI, or (b) Xtr — A but EPI is 
responsible for pairing by virtue of some peculiarities in 
the equations describing superconductivity - for instance 
because of non-Migdal vertex corrections. Related to 
this dilemma, it is worth mentioning that a similar puz- 
zling situation [Xtr ^ A) is realized in Ba^Ki-xBiOs 
(with Tc ~ 30 K), where optical measurements give 
Xtr ~ 0.1 — 0.3 28], while tunnelling measurements |29l | 
give A ~ 1. Note, in BoxKi-xBiOz there are no mag- 
netic fluctuations (or magnetic order) and no signs of 
strong electronic correlations. Therefore, EPI is favored 
as the pairing mechanism in Ba^Ki-xBiO^. It seems 
that in this compound long-range forces in conjunction 
with some nesting effects, may be responsible for this 
discrepancy. 

According to the above discussion, the microscopic the- 
ory of EPI and pairing in cuprates must be able to explain 
the following three important facts: (1) why Xtr ^ A; 
(2) why d-wave pairing is realized in the presence of 
strong EPI; (3) why is d-wave pairing robust in pres- 
ence of non-magnetic impurities and defects? Here we 
shall argue that the theory which is based on the exis- 
tence of the forward scattering peak (FSP) in EPI can 
account for both facts. We point out that d-wave pairing 
in cuprates is very robust in the presence of nonmagnetic 
impurities contrary to the prediction of the standard the- 
ory. We shall argue in what follows that this robustness 
is again due to strong correlations whic h g ive rise also 
to FSP in the impurity scattering 0, [ij. The FSP 
in EPI and in non-magnetic impurity scattering can be 
due to: (i) strong electronic correlations; (ii) the pres- 
ence of long-range Madelung potential in EPI, which is 
poorly screened for phonons vibrating along the c-axis; 
(iii) out of plain impurities and defects which are poorly 
screened. The concept of FSP in EPI and nonmagnetic 
impurity scattering potential due to strong correlations 
was for the first time put forward in ,30J, _31j, 32J. This 
approach, which was elaborated further by t he p resent 
author and his collaborators, see more in LZI, UJ, IS es- 
sential to explain ARPES in cuprates. 

This paper is organized as follows. In Sec. II we discuss 
direct and indirect experimental evidence for the impor- 
tance of the EPI in cuprates. The recent ARPES results 
are discussed in more detail. Section III is devoted to the 
systematic and controllable theory of strong correlations 
without using the slave boson (or fermion) technique. 
This theory treats strong correlations in terms of Hub- 
bard operators X"^ - we call it the X-method - which 
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describe the motion of the composite electrons. The lat- 
ter takes care that the double occupancy at a given lat- 
tice site is excluded. The central question in strongly 
correlated systems is - how to calculate the coherent part 
in the quasiparticle dynamics! This was achieved within 
the X-method by using a systematic and controllable 1/N 
expansion. In Sec. IV we discuss the renormalization 
of EPI by strong correlations. The 1/N expansion in- 
dicates that strong correlations lead to FSP in EPI. We 
stress that this result is straightforward in the X-method, 
while in the slave-boson method until now there have 
been no reliable calculations of EPI. We point out that 
the X-method predicts in the t-J model, that the EPI 
coupling becomes long-ranged in the presence of strong 
correlations. As the result of this specific renormaliza- 
tion, the EPI coupling constant in the d-channel reaches 
the value of that in the s-channel for some optimal dop- 
ing and below it. Furthermore, in conjunction with the 
residual Coulomb interaction this gives rise to d-wave 
pairing, since the residual Coulomb interaction is much 
stronger in the s-channel than in the d-channel. It is also 
shown that FSP appears in the nonmagnetic impurity 
scattering giving rise to robustness of d-wave pairing in 
cuprates. Section V applies the proposed theory to expla- 
nation of ARPES experiments, of the d-wave robustness 
in presence of impurities, of transport properties, etc. 
The problem of the pseudogap is discussed in Sec. VI 
from the point of view of the proposed theory. It turns 
out that the long-range EPI force gives rise to internal 
fluctuations of Cooper pairs, which together with phase 
fluctuations reduce strongly the mean-fleld critical tem- 
peratures. 

There are numerous experimental indications on the 
nanoscale inhomogeneities of HTSC oxides, for instance 
recent STM experiments [l^] . In that case the supercon- 
ducting order parameter is inhomogeneous, i.e. A(k, R) 
where R is the center of mass of Cooper pairs. In this 
review we shall not discuss this effects, which are in 
our opinion play secondary role in pairing mechanism. 
By understanding what is going on in homogeneous sys- 
tems we have a starting basis to study effects of inhomo- 
geneities. 



II. EXPERIMENTAL RESULTS IN FAVOR OF 
THE EPI 



The experimental situation which is related to the 
quasiparticle scattering in the normal and superconduct- 
ing state of cuprates is extensively analyzed in a number 
of reviews 0, |ll| |. Here we summarize briefly some ex- 
periments which are important for the microscopic theory 
of pairing in cuprates and which indicates importance 
of EPI in the quasiparticle scattering. Before this dis- 
cussion let us analyze briefly interplay of EPI and SFI 
in the pairing mechanism and their relative strength by 
taking into account three important experimental facts: 
(1) Tc in cuprates is high, i.e. T'^'''' ^ 160 AT; (2) 



the pairing in cuprates is d-wave like, i.e. A(k, k, 
As{k, oj) + Ad(w)(cos fcx — cos ky) with Ag < O-lAd, and 
(3) the EPI coupling constant is rather large, A « 1 — 2. 
Let us assume that d-wave pairing in cuprates is due to 
some non-phononic mechanism, for instance due to SFI. 
If EPI in cuprates would isotropic like in LTSC materi- 
als then it would be very detrimental for d-wave pairing 
as it is shown in by numerical studies of interplay 
of SFI and EPI by the Eliashberg theory. Here we shall 
elucidate this important problem in a physically plausi- 
ble and semi-quantitative way. In the case of dominating 
isotropic EPI in normal state, then near Tc the linearized 
Eliashberg equations have an approximative form 
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For pure d-wave pairing one has y(k, q, cj) = 
V{uj){coskj; — cosky){cosqx — cosqy) and A(k, cj) = 
A(a;)(cosfc;r — cosfcj,) what gives the equation for Tc 
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Here, T^o is the bare critical temperature due to the SFI 
mechanism and \& is the di-gamma function. At temper- 
atures near Tc one has T^p ~ 27rAepT and the solution 
of Eq. is 

Tc^Tcoe-^'-. 

The latter result means that when Aep = 1 — 2 and 
ymax _ ;^gQ ^ ^Yie bare Tco due to SFI must be very 
large, i.e. TcO = (400 - 1100) K, what is highly unreal- 
istic. However, this result, together with d-wave pairing 
in cuprates, imply that EPI in these materials must be 
strongly dependent on the momentum transfer, as it was 
flrst argued in ^^30]. Only in that case is EPI conform 
with d-wave pairing, either as its main cause or as a sup- 
porter of non-phononic mechanism such as SFI. In the 
following Sections we shall argue that the strongly mo- 
mentum dependent EPI is the main player in cuprates 
providing strength of the pairing mechanism, while the 
residual Coulomb interaction and SFI, although weaker, 
trigger it to d-wave pairing. 



A. Magnetic neutron scattering and spin 
fluctuation spectral function 

In a number of papers the experimental results for 
the pronounced imaginary part of the susceptibility 
/mx(k, fcz,a;), at the AF wave vector k = Q = (tTjTt), 
are interpreted as a support for SFI mechanism of pair- 
ing |l2j |. We briefly explain the inadequacy of such an 
interpretation. 
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(1) The breakthrough in understanding the role of 
spin-fluctuations came from magnetic neutron scattering 
experiments on YBa2Cu30e+x by Bourges group [33 |. 
They showed that the spectral function Imx''°'^'^^(k,uj) 
(the odd part in the bilayer system) is strongly depen- 
dent on hole doping. By varying doping there is a huge 
reconstruction of /mx'^°'^'')(Q, in the frequency inter- 
val which is important for superconducting pairing, say 
5 meV < u! < (60 — 70) meV. On the other hand, there 
is only a small variation of the T^. For instance, in un- 
derdoped YBa2Cu30e,g2 one has that Jtox^°'''^HQ) ^) 

(and S{Q) = /™d^Imx(°''^)(Q,^) - ^sfi {^)) is 
much larger than that in the near optimally doped 
YBa2Cu30Q,9T, i.e. 56.92(Q) > 5'6.97(Q), although the 
difference in their critical temperatures Tc is very small, 
i.e. Tc = 91 K in YBa2Cu30G,92 and = 92.5 K 
in YBa2Cu30Q,gj. The reader can convince himself di- 
rectly from the Fig^b). Such a large reconstruction of 
Imx''°''''HQ'^) but a negligible change in Tc in YBCO is 
a strong argument against the SFI mechanism of pairing 
in cuprates. Furthermore, the anti- correlation between 
the NMR spectral function /q = lim;j^n/mx(Q, ijj)/a; 
and Tc (where Tc ~ l//q, a > 0) see additionally 
disfavors SFI models of pairing i.e. the strength of 
pairing interaction is little affected by SFI. These results 
also mean that in spite of pronounced spin fluctuations 
around Q, the pairing quasiparticle self-energy is little 
affected by SFI. 

(2) The next very serious argument against SFI pair- 
ing mechanism is the smallness of coupling constant g^f . 
Namely, the real spin-fluctuation coupling constant ex- 
tracted from experiments is rather small g^f < 0.2 eV, in 

contrast to the large value {g[j^^'^'^ ~ 0.6 eV) assumed 
in the SFI theory by Pines and collaborators in order to 
get large Tc ~ 100 K. (Note that the pairing coupling in 
the SFI theory is As/ ~ g^/ ' ^^'^ realistic value of 

dsf < 0.2 eV it would produce Ag/ ~ 0.2 and very small 
Tc I K - see more in 0.) The large value for g[^^^^ 
gives resistivity much larger than what is observed 0. 
The upper limit for gsf{< 0.2 eV) is extracted: (i) from 
the width of the 41 meV magnetic resonance peak in the 
superconducting state |35l |: (ii) from the small magnetic 
moment {fi < 0.1 fj,^) found in the antiferromaenetic 
state of LSCO and YBCO. based on this fact in [Sg it 
was estimated As/ ^ 0.2; (iii) from ARPES experiments 
[21) that give As/ '--^ 0.3 - see also |3^ . 

(3) The SFI model is in some sense a phenomenolog- 
ical derivative of the t-J Hamiltonian. However, recently 
it was shown that there is no superconductivity in the t-J 
model at temperatures characteristic for cuprates [T^ . If 
it exists (at all!) Tc must be very low. 

(4) In the superconducting state the magnetic fluctu- 
ations are drastically changed, as expected for the sin- 
glet pairing state, since it induces a spin gap in mag- 
netic excitation spectrum of s-wave superconductors - 
see Fig.^c). However, the spectrum in superconducting 
state of cuprates is more complex due to d-wave pairing 
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FIG. 1: Magnetic spectral function Imx^~\^,i^)- (b) for 
YBa2Cu30e+x in the normal state at T = 100 K and at 
Q = (tt, 7r). 100 counts in the vertical scale corresponds to 
Xmax ~ SSO/il/eV; (c) for YBa2Cu30e+x in the supercon- 
ducting state a.t T — 5 K and at Q = (vr, vr). From Ref. 

a. 
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and due to some specificities of the band structure. For 
instance, at T < the sharp peak in /mx*^°'''*^(k, a;) 
is observed at ujreson = 41 meV and = {ir/a^n/a) 

- the resonance mode - of the fuUy oxygenated (opti- 
mally doped) Y Ba2Cu^OQ+x (x ~ 1, Tc w 92 K) - see 
the bottom of Fig. With increased doping, the res- 
onance mode becomes sharper and moves to higher fre- 
quencies (scahng with Tc), while its height decreases |^ 

- see Fig. ^c) . In optimally doped cuprates the spectral 
weight of the resonance peak is only 2-5 % of the total 
weight. Recently there were speculations in the literature 
that d-wave superconductivity in cuprates might be due 
to the resonance mode! This claim is untenable, since this 
mode can not cause superconductivity simply because its 
intensity near Tc is vanishingly small, thus not affecting 
Tc at all. If the magnetic resonance would be the origin 
for superconductivity (and high Tc) the phase transition 
at Tc would be of the first order, contrary to the experi- 
ments showing that it is of the second order. In fact, the 
resonance mode is a affected by superconductivity but not 
its cause. 

However, in spite of the fact that the strength of SFl 
is small, it can, together with other contributions of 
the residual Coulomb interaction, trigger d-wave pair- 
ing. However, the strength of pairing is due to EPI with 
FSP - see more below and in 0, [3- 



B. Dynamic conductivity and resistivity 

Optics measurements play very important role in 
studying the quasiparticle dynamics since this method 
probes bulk sample contrary to ARPES and tunnelling 
which probe tiny regions near the sample surface. 

A diversity of the results on the reflectivity measure- 
ments in cuprates were reported in the past. The dy- 
namic conductivity (t(w), which is not a directly mea- 
sured quantity but is derived from R{uj), was analyzed 
in terms of the quasiparticle self-energy E(fc,w) instead 
of the transport self-energy 'EtrikjUj). In fact 'E,{k,uj) ^ 
T,trik,U!) and this was the main reason for a number of 
erroneous conclusions. Here we discuss briefly the nor- 
mal state (7{u;) in the low frequency region to < 1 eV 
where the intra-band effects dominate the quasiparticle 
scattering. In the low ui regime, the processing of data in 
the metallic state of cuprates is usually done by using the 
generalized Drude formula for the in-plane conductivity 
ct(w) =ai(a;)+i(T2(cj) [H, [H 

uj\^ 1 
a^^{uJ) = ^^fT-^ — 7p{ ■ TTT • (3) 

An V tr\0J , 1 ) - iujmtr\'jJ) / moo 

i — a,b enumerates the plane axis, cop^a, Ttr{Lo,T) and 
mtr{i-o) are the plasma frequency, transport scattering 
rate and optical effective mass, respectively. For opti- 
mally doped HTSC systems the best fit for rtr{(jJ,T) is 
given by Ttrii^jT) w max{aT, (3lu} in a broad temper- 
ature (100 K < T < 2000 K) and frequency range (10 



meV < u! < 200 meV), where a,/? are of the order one. 
These results tell us that the quasiparticle liquid - which 
is responsible for transport properties in HTSC - is not a 
simple weakly interacting Fermi liquid. In that respect it 
is a well known fact that quasiparticles interacting with 
phonons at finite T are not described by the standard 
Landau-Fermi liquid, in particular at higher tempera- 
tures T > 0_d/5, since the scattering rate is of the order 
of quasiparticle energy, i.e. one has F ^ max(a;, T). Such 
a system is well described by the Migdal-Eliashberg the- 
ory whenever ujd Ep is fulfilled, which in fact treats 
quasiparticles beyond the original Landau quasiparticle 
concept. Note that even when the original Landau quasi- 
particle concept fails, the transport properties may be 
described by the Boltzmann equation, which provides a 
broader definition of the Landau-Fermi liquid. 

As we said, in a number of articles it was incorrectly 
assumed that r(cj,T) w rfr(w,T) holds in cuprates. 
Since the theory of EPI gives that T{uj,T) = const for 
Lu > ujph,ma.x (the maximal phonon frequency) and exper- 
iments give Ttr{uj, T) ^ uj for uj > tOph.ma.^, it was (incor- 
rectly) concluded that EPI in cuprates is small and can 
not be responsible for superconductivity and transport 
properties in the normal state. However, the thorough 
theoretical calculations of a{Lu) in a number of materi- 
als (including HTSC and LTSC systems) by including 
realistic EPI spectral function af^ ppF{uj) |3^], |3^, |40| | 
showed that as a rule one has: (1) Vtr{'-^,T) ^ V(uj,T) 
and (2) rtr(w,T) ~ w is linear in the broad region up 
to ujsat — -ImT.triiOsat) > Wp/i,max- Such a bchavior of 
rt,.(w,T) and T{u),T) is also realized in HTSC materi- 
als and in other metallic non-superconducting oxides; an 
extensive discussion is given in Q and Q . 

The main conclusion from the optics measurements is 
that EPI can explain the optical spectra of cuprates for 
Lo < 1 eV. Moreover, this analysis shows that the EPI 
coupling constant is rather large, i.e. 1 < Xepi ^ 2. 
The transport spectral function a^^{uj)F(uj) can be in 
principle extracted from the transport scattering rate 

Ttriuj,T = o) = — f dn{Lu -n)al{n)F{n) (4) 

^ Jo 

- see 0i 0- However, real measurements are performed 
at finite T(> Tc) where af^{u})F{u}) is the solution of 
the inhomogeneous Fredholm integral equation of the 
first kind. This inverse problem at finite t emp eratures 
in cuprates was studied in [s^ (see also UJ), where 
the smeared structure of a^j.{uj)F{uj) in Y Ba2Cu^O-j-x 
was obtained, in a qualitative agreement with the shape 
of the phonon density of states F{uj). At finite T the 
problem is more complex because the fine structure of 
q;j^((jj)F((jj) gets blurred, as the calculations in 41] con- 
firm. The latter showed that a(^(a;)F(w) ends up at 
^max ~ 70 — 80 meV , which is the maximal phonon fre- 
quency in cuprates. This result indicates strongly that 
the EPI in cuprates is dominant in the IR optics. Strictly 
speaking, this is still not a definitive proof for the EPI 
since the solution of the Fredholm equation of the first 
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kind is extremely sensitive to the input, because the un- 
known function {a'jj,{Lij)F(oj) which is a positive quan- 
tity) is under integral and smoothed by some kernel. It 
is known that this smoothing can miss some fine details 
of the solution. 

We point out that if R{uj) (and are due to some 

other bosonic process with large frequency cutoff lOc in 
the spectrum, as it is the case with SFI scattering where 
aj^((jj)F((jj) ~ gl^ImXsi'-jj) and ujc ~ 400 meV , the ex- 
tracted a1j.{u})F{Lo) should end up at this high lUc- This 
was done in 42] for T > Tc, where it was found that 
afj.{uj)F(uj) reaches negative values in a broad energy re- 
gion 70 meV < tu < coc- These results rule out SFI as 
a dominant scattering mechanism in cuprates and favor 
EPI. Recent 

We stress that the extraction of Ttr from -R(w) is a 
rather subtle procedure, since it depends on the assumed 
value of dielectric constant Eoo- For instance, if one takes 
£oo = 1 then is linear up to very high cu, while for 
£oo > 1 the linearity of saturates at lower oj. For 
instance, the extracted r^^(a;,T) in is linear up to 
very high uj. Since in these papers there is no informa- 
tion (!) on Eoo, the too linear behavior might mean that 
the ion background and interband transitions (contained 
in Eoc) are not properly taken into account in these pa- 
pers. We stress again, that the behavior of Ttrii^j) is 
linear up to much higher frequencies for Eoo = 1 than for 
£oo ~ 4 — 5 - the characteristic value for HTSC, giving 
a lot of room for inadequate interpretations of results. 
Note that recent elipsometric measurements on YBCO 
[4^ confirm the earlier results |4^ that £00 > 4 and that 
F^^ saturates at lower frequency than it was the case in 
Ref. ji^l . We stress again that a reliable estimate of the 
value and of the cu, T dependence of Ttrii^j) and rn(uj) can 
be done, not from the reflectivity measurements [43j, but 
from elipsometric ones only (4^] , 0| . 

In concluding this part we stress two facts: (1) The 
large difference in the uj,T behavior of Tf/'^ {u!,T) and 
r^^^ {uj,T) is not a specificity of cuprates but it is real- 
ized also in a number of LTSC materials. In fact this 
is a common behavior even in simple metals, such as Al 
or Pb where T^^' {u!,T) saturates at much lower (De- 
bay) frequency than F^^^(u;,r) and F*;^^^(w,T) do - 
see more in Q, Q and references therein. In that re- 
spect the difference between simple metals and cuprates 
is in the scale of phonon frequencies, i.e. i^ai&x ^ 100 
K in simple metals, while oj^g^^ ~ 1000 K in cuprates. 
Having in mind these well established and well under- 
stood facts, it is very surprising that even nowadays, 
twenty years after the discovery of cuprates, the essen- 
tial and quantitative difference between F and Fi^ is ne- 
glected in the analysis of experimental data. For in- 
stance, by neglecting the pronounced (qualitative and 
quantitative) difference between Ttr{uJ,T) and T{uj,T), 
in the recent papers |43l | were made far reaching, but 
unjustified, conclusions that the magnetic pairing mech- 
anism prevails; (2) It is worth mentioning that quite sim- 
ilar (to cuprates) properties, of cr(w), ^(w) and p(T) were 



observed in experiments |46| on isotropic metallic oxides 
Lao.^Sro^^CoOs and Cao.^Sro ^RuO^. We stress that in 
these compounds there are no signs of antiferromagnetic 
fluctuations (which are present in cuprates) and the pe- 
culiar behavior is probably due to the EPI. 

In that respect, recent experiments related to the opti- 
cal sum-rule with the transfer of the spectral weight from 
high to low frequency is an additional example of contro- 
versies in this field coming from inadequate interpreta- 
tion of experimental results. We stress that the theory 
based on EPI can explain the transfer of spectral weight 
by taking into account strong the co and T dependence 
of F^(u;,r) due to the EPI - for details see [H, Q 
and 113 . 

The temperature behavior of in-plane resistivity 
Pa-bi'^) is a direct consequence of quasi-2D motion of 
quasiparticles and of inelastic scattering which they ex- 
perience. At present, there is no consensus on the origin 
of linear temperature dependence of the in-plane resistiv- 
ity Pa_fc(T) in the normal state. Many researchers (er- 
roneously) believe that such a behavior can not be due 
to the EPI. The inadequacy of this claim was already 
demonstrated by analyzing the dynamic conductivity 
(t{uj). It is well-known that at temperatures T > 0d/5 
and for the Debay spectrum one has 

p(T)^8.2Af/^=p'r. (5) 

In cuprates the reach and broad spectrum of a^^{u))F{ul) 
favors such a linear behavior in a broader T region. 
The measured transport coupling constant Xtr contains 
in principle all scattering mechanisms, although usually 
some of them dominate. For instance, the proponents of 
SFI mechanism assume that Xtr is entirely due to scatter- 
ing on spin fluctuations. However, by taking into account 
specificities of cuprates the experimental results for the 
in-plane resistivity Pa^ti'^) *^an be satisfactory explained 
by the EPI mechanism. From tunnelling experiments 
one concludes that A « 2 — 3 and if one assumes that 
Xtr ~ X and ujpi > (3 — 4) eV (the value obtained from 
the band-structure calculations) then Eq. ||SJ) describes 
the experimental situation rather well. The plasma fre- 
quency u>pi which enters Eq.ljSj can be extracted from 
optic measurements (ujpi^ex), i-e. from the width of 
the Drude peak at small frequencies. However, since 
At,, w 0.25ujli{eV)p'{nncm/K) there is an experimental 
constraint on Af^- This gives 46j that uipi « (2 — 2.5) eV 
and p' w 0.6 in the oriented YBCO films, while p' « 0.3 
in single crystals of BSCCO. These results impose a limit 
on Xtr ~ 0.4 - 0.8. 

The above analysis implies, that in order to explain 
p(T) with small Xtr and high Tc (which needs large A) 
by EPI it is necessary to have Xtr < (-^/S)- This means 
that in cuprates EPI is, due to some reasons, reduced in 
transport properties where Xtr ^ A. This reduction of 
and At,- means that the A and At,- contain renormal- 
ization (with respect to the LDA results) due to various 
quasiparticle scattering processes and interactions, which 
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do not enter in the LDA theory. In subsequent chapters 
we shall argue that the strong suppression of Xfj may 
have its origin in strong electronic correlations i j 

M 

In conclusion, optics and resistivity measurements in 
normal state of cuprates are much more in favor of EPI 
than against it. However, some intriguing questions still 
remains to be answered: (i) what are the values of Xtr 
and LOpi; (ii) what is the reason that Xtr «C A is realized 
in cuprates; (iii) what is the role of Coulomb scattering 
in a{uj) and p{T). Later on we shall argue that ARPES 
measurements in cuprates give evidence for a contribu- 
tion of Coulomb scattering at higher frequencies, where 
r(w) « Fo + AcW for uj > LjP^^^ with Ac « 0.4. So, 
in spite of the fact that EPI is suppressed in transport 
properties it is sufficiently strong that it dominates in the 
self-energy in some frequency and temperature range. It 
is possible that at higher temperatures Coulomb scatter- 
ing dominates p{T), which certainly does not disqualify 
EPI as the pairing mechanism in cuprates. For better 
understanding of p{T) we need a controllable theory for 
Coulomb scattering in strongly correlated systems, which 
is at present lacking. 



C. Electronic and phonon Raman scattering in 
cuprates 

(a) Electronic Raman scattering in various cuprates 
show a remarkable correlation between the Raman cross- 
section S'exp(w) and the optical conductivity aa-b(y), i-e. 
Sexpl'^) ~ o'a-b(w) 0- Previously it was demonstrated 
that (Ta_b(w) depends on the transport scattering rate 
Ft,.(w, T). We have also demonstrated that the EPI with 
the very broad spectral function a'^F{u!) (see Fig. 2 be- 
low) explains in a natural way the w, T dependence of 
(Ta_b(cij) and Ttr{^,T). Therefore, the Raman spectra 
in cuprates can be explained by the EPI in conjunc- 
tion with strong correlations. This conclusion is sup- 
ported by calculations of the Raman cross-section [47| 
which take into account the EPI with a^FioS) extracted 
from the tunnelling measurements on Y Ba2CuzOQ+x 
and Bi2Sr2CaCu20s+x E3; 0- Quite similar (to 
cuprates) properties of the electronic Raman scattering 
(besides a{uj), R{ijj) and p{T)) were observed in exper- 
iments |46| on isotropic metallic oxides Lao.sS'ro.sCoOs 
and Cao.^Sro.^RuOs where there are no signs of antifer- 
romagnetic fluctuations. 

(b) Two important experimental results related to the 
phonon Raman scattering give evidence for strong EPI 
in cuprates: (i) there is a pronounced asymmetric line- 
shape (of the Fano resonance) in the metallic state. For 
instance in Y Ba2Cu^0'! two Raman modes at 115 cm^^ 
(Ba dominated mode) and at 340 cm^^ (O dominated 
mode in the Cu02 planes) show pronounced asymmetry 
which is absent in Y Ba2Cu^O(,. This result points to the 
strong interaction of Raman active phonons with contin- 
uum states (quasiparticles) ; (ii) The phonon frequencies 



for some Aig and Big are strongly renormalized in the 
superconductingstate, pointing again to the large EPI - 
see more in 0,^3- 



D. Tunnelling spectroscopy and Eliashberg 
spectral function 



By measuring the I-V characteristic, i.e., the tunnelling 
conductance G{V) = dl/dV, one can determine the 
Eliashberg spectral function a'^F{uj) which determines 
the pairing coupling constant X — 2 J dLua'^F{uj)/uj. A 
number of experiments 17] on tunnelling- and break- 
junctions in cuprates gave a^F(w) which is spread over a 
broad range of frequencies (0 < w < 80 meV) and whose 
maxima coincide with the maxima in the phonon density 
of states F{u!) - see Fig. 2. This is a rather strong proof 
of importance of the EPI in pairing. Moreover, most of 
these experiments give the range of the coupling constant 
1 < A < 2.5. 

Since practically all phonon modes contribute to A - 
see Fig. 2, this means that on the average each particu- 
lar phonon mode is moderately coupled to electrons thus 
keeping the lattice stable. Additionally, they have found 
that some low-frequency phonon modes corresponding to 
Cu, Sr and Ca vibrations are rather strongly coupled to 
electrons. The high frequency oxygen vibrations along 
the c-axis interact with quasiparticles appreciably too. 
These results confirm the importance of axial phonon 
modes in which the change of the Madelung energy is 
involved, thus supporting the idea conveyed through this 
article of the importance of the ionic Madelung energy in 
the EPI interaction of cuprates. 

Recent STM experiments ^3 confirmed partly the re- 
sults of tunnelling measurements. They find near the 
anti-nodal point strong EPI with 40 meV phonons and 
from Fourier-transformed spectra the exclude the mag- 
netic resonance as a possible origin. 

In conclusion, the common results for tunnelling mea- 
surements in cuprates, including Bai-xK^BiO^ too [48l| . 
is that no particular mode can be singled out in the 
spectral function a'^F{uj) as being the only one which 
dominates in pairing mechanism. This important result 
means that the high Tc is not attributable to a particu- 
lar phonon mode in the EPI mechanism, since all phonon 
modes contribute to A. Having in mind that the phonon 
spectrum in cuprates is very broad (up to 80 meV ), 
then the large EPI coupling constant (A « 2) in cuprates 
is not surprising at all. We stress, that compared to neu- 
tron scattering experiments the tunnelling experiments 
are superior in determining the EPI spectral function 
a'^F{ui) - see more in 0. 
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FIG. 2: The spectral function cPf{uj) obtained from mea- 
surements of G{y) by various groups on various junctions. 
The generahzed density of states GPDS for Bi2212 is plotted 
at the bottom. From Ref. \vi\. 



E. Isotope effect for various doping 

It is well known that in the pure EPI pairing mecha- 
nism the total isotope coefficient a is given by 



E 



dhiTr 



(6) 



where M^^^ is the mass of the i-th element in the p- 
th crystallographic position. From this formula one can 
deduce that the relative change of Tc, 5Tc/Tc, for heav- 
ier elements is rather small - for instance it is 0.02 for 
135^^ ^138 0.03 for 63Cy ^65 (j^ ^^^^ g.O? for 

^^^La ^139 This means that measurements of ai for 
heavier elements are at/or beyond the ability of exper- 
imental techniques. Therefore most isotope effect mea- 
surements were done by substituting light atoms ^^O by 
^^O. In that respect one should have in mind the tun- 
nelling experiments discussed above, which tell us that 



0.1 — 0.2 in Lai,s5Sro,i5Cu04 with 
w 0.03 - 0.05 in Bi2Sr2CaCu208 



practically all phonons contribute to the Eliashberg pair- 
ing function a'^F{'k, uj), and that oxygen modes give only 
moderate contribution to T^,. So, small oxygen isotope 
effect does not exclude the EPI mechanism of pairing. 

Measurements of the isotope coefficients ao (and acu) 
were performed on various hole-doped and electron- 
doped cuprates. A brief summary of the main results is 
as follows HI: (1) The O isotope coefficient ao strongly 
depends on the hole concentration in the hole-doped ma- 
terials; in each family of cuprates {Y Ba2Cu30j^x, or 
La2-xSrxCuOi etc.) a small oxygen isotope effect is 
observed in the optimally doped (maximal Tc) samples. 
For instance ao ~ 0.02 — 0.05 in Y Ba2Cu^Oj with 
Tcmax ~ 91 K, ao ~ 

.max 

~ 35 K; ao 

with Tcmax ~ 76 K; ao ~ 0.03 and even negative 
(-0.013) in BiaSrzCaaCusOio with Te^^ax « HO K; the 
experiments on Tl2Can-iBaCun02n+4 [n = 2, 3) with 
Tc,max ~ 121 K are still unreliable and ao is unknown; 
ao < 0.05 in the electron-doped {Ndi-xCex)2Cu04 
with Tcmax ~ K . (2) For hole concentrations away 
from the optimal one, Tc decreases while ao increases 
and in some cases reaches large value ao ~ 0.5. This 
holds not only for parent compounds but also for systems 
with substitutions, like {Yi-x-y^T^x Cay)Ba2Cu30j, 
Yi^yCayBa2Cu404 and Bi2Sr2Cai-xYxCu20s- Note 
that the decrease of Tc is not a prerequisite for the in- 
crease in ao- This became clear from the Cu substituted 
experiments Y Ba2{Cui-xZnx)30'/ where the decrease of 
Tc (by increasing the Zn concentration) is followed by 
only a small increase of ao- Only in the case of very 
low Tc < 20 K, ao becomes large, i.e. ao > 0.1. (3) 
The largest ao is obtained even in the optimally doped 
compounds like in systems with substitution, such as 
Lai,s5Sro,i5Cui-xMx04, M = Fe,Co, where ao ~ 1.3 
for X « 0.4 %. (4) In La2-xMxCu04 there is a Cu 
isotope effect which is of the order of the oxygen one, 
i.e. acu ~ cto giving acu + ao ~ 0.25 — 0.35 for op- 
timally doped systems {x — 0.15). In the case when 
X — 0.125 with Tc <C Tcmax one has acu ~ 0.8 — 1 with 
acu + OLo « 1.8. The appreciate copper isotope effect in 
La2-xMxCu04 tells us that vibrations of other than oxy- 
gen ions could be important in giving high Tc. The latter 
property is more obvious from tunnelling measurements 
which are discussed above. (5) There is a negative Cu iso- 
tope effect in the oxygen-deficient system Y Ba2Cu30i-x 
where acu is between —0.14 and —0.34 if Tc lies in the 
60 K plateau. (6) There are reports on a small negative 
ao in some systems like Y Sr2Cu30j with ao ~ —0.02 
and in BSCCO -222?, {Tc = 110 K) where ao ~ -0.013 
etc. However, the systems with negative ao present con- 
siderable experimental difficulties, as it is pointed out in 

The above enumerated results, despite experimental 
difficulties, tell us that the EPI interaction is involved 
in the pairing mechanism of cuprates and that the rela- 
tion of Tc to ionic masses is highly nontrivial and non- 
BCS-like. By assuming that the experimental results on 



10 



the isotope effect reflect an intrinsic property of cuprates 
one can rise a question: which theory can explain these 
results? Since at present there is no consensus on the 
pairing mechanism in HTSC materials there is also no 
definite theory for the isotope effect. Besides the calcu- 
lation of the coupling constant A any microscopic the- 
ory of pairing is confronted also with the following ques- 
tions: (a) why is the isotope effect small in optimally 
doped systems and (b) why a increases rapidly by fur- 
ther under(over)doping of the system? The answer to 
this question is additionally complicated by the quasi-2D 
structure of cuprates and the closeness of these materi- 
als to the Mott-Hubbard isolator. It is known that the 
later gives rise to the smallness of the phase stiffness 
of the superconducting order parameter A =| A | e*'''. 
In 2D superconductors the energy increase due to phase 
fluctuations in the conduction plane is given by [49j 

Fphase{V^,T) = ^K^J d?x{V^f, (7) 

where — h^ns{T)/Am* is the 2D stiffness, Ug is the 
2D superconducting charge density and m* « 2me is the 
effective mass. In the case of the Berezinskii-Kosterlitz- 
Thouless transition (vortex-antivortex unbinding) one 
has Tc = {-K /2)Ks{Tc), where Ks is the renormalized 
stiffness. In that case the isotope effect is completely 
different from the prediction of the ECS theory. In that 
respect the very careful experiments on isotope effect of 
the penetration depth in the a-b plane of YBCO show 
a significant increase upon the substitution O^^ O^*, 

i.e. (AAab/ Aafc) = {^^Xab "^^ Xab)l^''Xab = 2.8 % at 4 

K jl3|. Since Aafc ~ nn* this could mean that there is a 
nonadiabatic increase of the effective mass to* . 



others in ARPES experiments. After a drastic improve- 
ment of sample quality in other families of HTSC materi- 
als, the ARPES technique has became a central method 
in theoretical considerations. Potentially, it gives valu- 
able information on the quasiparticle Green's function, 
i.e. on the quasiparticle spectrum and life-time effects. 
The ARPES can indirectly give information on the mo- 
mentum and energy dependence of the pairing poten- 
tial. Furthermore, the electronic spectrum of the (above 
mentioned) cuprates is highly quasi-2D which allows an 
unambiguous determination of the initial state momen- 
tum from the measured final state momentum, since the 
component parallel to the surface is conserved in pho- 
toemission. In this case the ARPES probes (under some 
favorable conditions) directly the single particle spectral 
function A{\i,Lu). In the following we discuss only those 
ARPES experiments which give evidence for the impor- 
tance of the EPI in cuprates - see more in [21| . 

The photoemission measures a nonlinear response 
function of the electron system, and under some con- 
ditions it is analyzed in the so-called three-step model, 
where the total photoemission intensity /tot(k, w) « I • 
I2 ■ I3 is the product of three independent terms: (1) / 
- describes optical excitation of the electron in the bulk; 
(2) I2 - the scattering probability of the travelling elec- 
trons; (2) I3 - the transmission probability through the 
surface potential barrier. The central quantity in the 
three-step model is /(k, uj) and it turns out that it can 
be written in the form for k = k|| 21] 

I{k,u;)c,Io{k,v)f{u;)Aik,u;) 
/o(k,w)^|(^. |pA IV,) P 



F. Angle-resolved photoemission spectroscopy in 
cuprates 

ARPES is nowadays a leading spectroscopy method 
in the solid state physics The method consists in 

shining light (photons) with energies between 20 — 1000 
eV on the sample and by detecting momentum (k)- and 
energy(a;)-distribution of the outgoing electrons. The 
resolution of ARPES has been significantly increased in 
the last decade with the energy resolution of AE « 2 
meV (for photon energies ~ 20 eV) and angular reso- 
lution of Ae w 0.2°. The ARPES method is surface 
sensitive technique, since the average escape depth (lesc) 
of the outgoing electrons is of the order of l^sc 10 A, 
of course depending on the energy of incoming photons. 
Therefore, very good surfaces are needed in order that 
the results be representative of bulk samples. In that re- 
spect the most reliable studies were done on the bilayer 
Bi2Sr2CaCu20s {Bi2212) and its single layer counter- 
part Bi2Sr2CuOQ {Bi22Ql), since these materials con- 
tain weakly coupled BiO planes with the longest inter- 
plane separation in the cuprates. This results in a nat- 
ural cleavage plane making these materials superior to 



A{k, uj) = — ' M2 I v2n V 

TT - ^(k) - El (k, uj)Y + ^2 (k, w) 

Here, {tpj \ pA 1-0,) is the dipole matrix element which 
depends on k, polarization and energy v of the incoming 
photons. The knowledge of the matrix element is of a 
great importance and its calculation from first principles 
was done carefully in f{uj) = 1/(1 + exp{a;/T}) is 
the Fermi function, A{k,uj) = —lmG{k,uj)/TT, G{k,uj) 
and I](k, w) = T^i{k,uj) + iE2(k,w) are the spectral 
function, the quasiparticle Green's function and the self- 
energy, respectively. 

In some period of the HTSC era there were a number 
of controversial ARPES results and interpretations, due 
to bad samples and to the euphoria with exotic theories. 
For instance, a number of (now well) established results 
were questioned in the first ARPES measurements, such 
as: the shape of the Fermi surface, which is correctly pre- 
dicted by the LDA band-structure calculations; bilayer 
splitting in Bi2212, etc. 

We summarize here some important ARPES results 
which were obtained recently, first in the normal state 
[2ll |: (Iat) There is a well defined Fermi surface in the 
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metallic state - with the topology predicted by the LDA 
but the bands are narrower than the LDA ones; (2^?) the 
spectral lines are broad with | S2(k, w) \^ uj (or ^ T for 
T > uj); (Sat) there is a bilayer band splitting in Bi2212 
(at least in the over-doped state); {4^) at temperatures 
Tc < T < T* and in the under-doped cuprates there is a 
d-wave like pseudo-gap Apg(k) Apg^Q^cos — coaky) 
in the quasiparticle spectrum; (Sat) the pseudo-gap Ap^^o 
increases by lowering doping; {6n) the ARPES self- 
energy gives clear evidence that the EPI interaction, with 
the characteristic phonon energy Wp/i, is rather strong. 
At T > Tc there are kinks in the quasiparticle dispersion 
"^(Ck) ill the nodal direction (along the (0, 0) — (tt, tt) line) 
at w^"^"^ ~ (60 - 70) meV 3], see Fig. El and near the 
anti-nodal point (tt, 0) at 40 meV - see Fig. 13 {7n) 
The holes couple practically to the whole spectrum of 
phonons. For instance, at least three group of phonons 
(interacting with holes) were extracted recently from the 
ARPES self-energy in La^-xSr^CuOA 113 ■ (8^) The 
EPI coupling constant which is obtained by comparing 
the slope duj/d^^^ of the quasiparticle energy uj(S.k) at 
very small | Ck l< ^ph and | |> ujph gives Xept « 1-2, 
while the Coulomb part is Ac « 0.4. {9n) Recent ARPES 
spectra in the optimally doped Bi2212 near the nodal 
and anti-nodal point 24J show a pronounced isotope ef- 
fect in the real part of the quasiparticle self-energy, thus 
pointing out the important role of the EPI. 

The corresponding ARPES results in the supercon- 
ducting state are the following |23|: (Is) there is 
an anisotropic superconducting gap in most HTSC 
compounds, which is predominately d-wave like, i.e. 
Asc(k) ~ Ao(cosA:a; — cosky) with 2Ao/Tc ~ 5 — 6; 
(2s) there are dramatic changes in the spectral shapes 
near the anti-nodal point (tt, 0), i.e. a sharp quasiparticle 
peak develops at the lowest binding energy followed by 
a dip and a broader hump, giving rise to the so called 
peak- dip-hump structure; (3s) the kink in the quasipar- 
ticle energy around the nodal-point and at (60 — 70) 
meV is surprisingly not-shifted in the superconducting 
state. To remind the reader, in the standard Eliashberg 
theory the kink should be shifted to LOph + Aq at any 
point at the Fermi surface. (4s) the anti-nodal kink at 

by Ao, 

Aq = (25 — 30)TOeI^. In the following sections we shall 
argue that all these properties can be explained by EPI 
with FSP and the coupling constant A > 1. 



40 meV is shifted in the superconducting state 
Ao = (65 — 70)meV since 
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FIG. 3; (Top) Quasiparticle dispersion of Bi22l2, Bi220l 
and LSCO along the nodal direction, plotted vs the momen- 
tum k for (a) — (c) different doings, and (d) — (e) different T; 
black arrows indicate the kink energy; the red arrow indicates 
the energy of the q = (tt, 0) oxygen stretching phonon mode; 
inset of (e)- T-dependent E' for optimally doped Bi2212; (/) 
- doping dependence of A' alo ng ( 0, 0) — (tt, tt) for the differ- 
ent HTSC oxides. From Ref. [20|l . (Bottom) Quasiparticle 
dispersion E{k) in the normal state (al, bl, c), at 107 K and 
115 K, along various directions (p around the anti-nodal point. 
The kink at E = 4QmeV is shown by the horizontal arrow. 
(a2 and b2) is E{k) in the superconducting state at 10 K with 
the shifted kink to 70meV . (d) kink positions as a function 
of (j) in the anti- nodal region. From Ref. |2^. 



III. THEORY OF STRONG CORRELATIONS 
WITHOUT SLAVE BOSONS 

The well established fact is that strong electronic cor- 
relations are pronounced in cuprates. However, the LDA 
theory fails to capture effects of strong correlations by 
treating them as a local perturbation only. This is an un- 
realistic approximation for cuprates, where strong corre- 



lations introduce non- locality in charge interactions. One 
of the shortcomings of the LDA is that in the half-filled 
case (with n — 1 and one particle per lattice site) it pre- 
dicts metallic state only, thus missing the existence of the 
Mott insulating state. In the latter, particles are localized 
at lattice sites independently of the (non)existence of the 
AF. This localization is due to the large Coulomb repul- 



12 



sion U at given lattice site, i.e. U ^ W where W is the 
band width. Some properties in the metaUic state can not 
be described by the simple canonical Landau-Fermi liq- 
uid concept. For instance, recent ARPES photoemission 
measurements on the hole-doped samples show a well de- 
fined Fermi surface in the one-particle energy spectrum, 
which contains \ — 5 electrons in the Fermi volume {5 is 
the hole concentration), but the band width is (2 — 3) 
times smaller than the LDA band structure calculations 
predict |23| • The " 1 — (5" behavior of the Fermi volume 
is consistent with the Luttinger theorem. However, ex- 
perimental data on the dynamic conductivity (spectral 
weight of the Drude peak). Hall measurements etc. in- 
dicate that in transport properties a low density of hole- 
like charge carriers (which is proportional to 5) partic- 
ipates predominantly. These carriers experience strong 
scattering and their inverse lifetime is proportional to 
the temperature (at T > Tc) as we discussed earlier. It 
is worth mentioning here that the local moments on the 
Cu sites, which are localized in the parent AF compound, 
are counted as a part of the Fermi surface area when the 
system is doped by small concentration of holes in the 
metallic state. The latter fact gives rise to a large Fermi 
surface which scales with the number (per site) of elec- 
trons 1 — (5. At the same time the conductivity sum-rule 
is proportional to the number of doped holes 5, at least 
in underdoped systems, i.e., /p"" dijj(Ji{w) ~ 5) instead of 
1 — 5 as in the canonical Landau- Fermi liquid. These two 
properties tell us that in cuprates we deal with a cor- 
related state, and the latter must be due to the specific 
electronic structure of cuprates. The common ingredient 
of all cuprates is the presence of Cu atoms. In order 
to account for the absence of Cm^+ ionic configuration 
(the charge transfer Cm^+ Cu^^ costs large energy 
?7 ^ 10 eV, i.e. the occupation of the Cu site with two 
holes with opposite spins is unfavorable) P. W. Ander- 
son proposed in his remarkable paper |5ll | the Hubbard 
model as the basic model for quasiparticle properties in 
these compounds. For some parameter values it can be 
derived from the (minimal) microscopic three-hand model 
in some parameter regime. Besides the hopping tpd be- 
tween the d-orbital of Cu and p-orbital of O ions (as well 
as tpp between O ions) - the Emery model, it contains 
also the strong Coulomb interaction Ucu on the Cu ions 
as well as interaction between p- and d-electrons. The 
main two parameters are Ucu ~ (6 — 10) eV and the 
charge transfer energy Ap^ = — Cp ^ (2.5 — 4) eV , 
where e^, are bare energies of the d- and p levels, re- 
spectively. In cuprates the case Ucu >> ^pd is realized, 
i.e. they belong to the class of charge transfer mate- 
rials. This allows us to project the complicated three- 
band Hamiltonian onto the low-energy sector, and to ob- 
tain an effective single-band Hubbard Hamiltonian with 
an effective hopping parameter t and the effective repul- 
sion Ueff ~ Aprf. It turns out that the case Ueff » t 
is realized, if \d > t = t'^JApd- In fact in HTSC 
cuprates the latter is not quite realized and probably 
the slightly weaker condition Ap^ > t is more appro- 



priate, since tpd « 1.3 eV and Ap^ « (2 — 3) eV. Nev- 
ertheless, there is a widespread believe that the suitable 
effective and minimal Hamiltonian, i.e. the toy model 
for low- energy physics of cuprates, is based on the Hub- 
bard model which additionally comprises the long-range 
Coulomb interaction Vc and EPI Vepi - see 01 

H ^ - Y^ tijcl^Cj„ + Ueff Ui^n^i + Vc + Vepi- (9) 

Some experiments imply that the effective repulsion 
should be of the order of Ueff « 4 eV, while the near- 
est neighbor and next-nearest neighbor hopping t and 
t' , respectively are estimated to be t = 0.3 — 0.5 eV 
and t'/t equal -0.15 in LSCO and -0.45 in YBCO. 
Since (C/g/ f/t) ^ 1 the above Hamiltonian is again in the 
regime of strong electronic correlations, where the double 
occupancy of a given lattice site is strongly suppressed, 
i.e. {nii^riii) <^ 1. The latter restricts charge fluctuations 
of electrons (holes) on a given lattice site, since rii = 0, 1 
is allowed only, while processes with rii — 2 are (practi- 
cally) forbidden. Note, that in (standard) weakly corre- 
lated metals all charge fluctuation processes (n^ = 0, 1, 2) 
are allowed, since U <^W in these systems. Due to sup- 
pression of double occupancy one expects that the screen- 
ing properties in cuprates are strongly affected - this will 
be demonstrated below. From the Hamiltonian in Eq. 
, which is the 2D model for the low-energy physics in 
the Cu02 plane, comes out that in the undoped system 
there is one particle per lattice - the so called half-filled 
case (in the band language) with {ui) = 1. It is an in- 
sulator because of large U and even antiferromagnetic 
insulator at T = The effective exchange interaction 
(with J = At^/U) between spins is Heisenberg-like. By 
doping the system with holes (the hole concentration is 
i5(< 1)) means that particles are taken out from the sys- 
tem in which case there is on the average {rii) = 1 — 5 
particles per lattice site. Above some (small) critical dop- 
ing 5c (which is of the order 0.01 in doped cuprates) the 
AF order vanishes and the system becomes a strongly 
correlated metal. For some optimal doping 5op{'^ 0.1) 
the system is metallic with a large Fermi surface and can 
exhibit even high-Tc superconductivity in the presence 
of the EPI, as it will demonstrated below. The latter in- 
teraction and its interplay with strong correlations is the 
central subject in the following sections. 

A. Strong correlations in the t-J model in terms of 
Hubbard operators X"^ 

Since Ueff >> t one can let with a good accuracy 
Ueff — > oo, i.e. the system is in the strongly corre- 
lated regime where the doubly occupancy = 2 is ex- 
cluded. One of the ways to cope with such strong cor- 
relations is to introduce the Hubbard operators X"^ 
which automatically take into account the exclusion of 
double occupancy. The creation and annihilation Hub- 
bard operators and X^" = (XfO)t, with XfO = 



13 



c|g.(l — Tii^-a), are fermion-like and they fulfill the condi- 
tion that rii^a- + ni _^ < 1 on each lattice site. The latter 
means that there is no more than one electron (hole) 
per lattice site, i.e. the double occupancy is forbidden. 
The bosonic-like Hubbard operators X^^"^ = i°X°''^ 
(with (Ti ^02) describe spin fluctuations at the i-th site. 
Here, the spin projection parameter a =t, i ( —o =i, t) 
and the operator X"" has the meaning of the electron 
(hole) number on the i-th site. In the following we shall 
use the convention that when Xf"' | 1) = 1 | 1) there 
is a fermionic particle (electron) on the j-th site, while 
for Xf"" I 0) = I 0) the site is empty, i.e. there is 
a hole on it. It is useful to introduce the hole number 
operator X'f = X^^X^^, i.e. if X°° | 0) = 1 | 0) the 
j-th site is empty, i.e. there is one hole on it, while for 
X™ I 1) = I 1) it is occupied by an electron and there 
is no hole. These operators fulfill the non-canonical com- 
mutation relations 



= s.. 



7/3 



(10) 



Here, a,/?, 7, A — 0,a and Sij is the Kronecker sym- 
bol. The (anti)commutation relations in Eq. H10() 
are rather different from the canonical Fermi and Bose 
(anti)commutation relations. Since Ueff = 00 the double 
occupancy is excluded, i.e. hiifiii \ ip){= XP |ti)) = 0, 
and by construction the X"'' operators satisfy the local 
constraint (the completeness relation) 
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(11) 



This condition tells us that at a given lattice site there is 
either one hole (X™ | hole) — 1 | hole)) or one electron 
(Xf" I elec) = 1 I elec)). Note, if Eq. ^ is obeyed, 
then both commutation and anticommutation relations 
hold also in Eq. (|l()|l at the same lattice site, which is 
due to the projection properties of Hubbard operators 

For further purposes, i.e. for studying low-energy ex- 
citations in a controllable way, the number of spin pro- 
jections is generalized to be X , i.e. a = 1,2,...X. By 
projecting out the double occupied (high energy) states 
from the Hamiltonian in Eq. © one obtains the gener- 
alized t-J model [111 



of almost-localized electrons. Hz contains three-site 
hopping which is believed not to be important and is 
usually omitted. For effects related to charge fluctua- 
tion processes it is plausible to omit it, while for spin- 
fluctuation processes this approximation may be ques- 
tionable. The spin and number operators are given by 
S — X^^^{a)ai92^i'^^'j where summation over 

bar indices is assumed [flj . 

The rather "awkward" commutation relations of X"^^ 
were one of the reason that most researchers in the field 
prefer to use some representations of the Hubbard op- 
erators X"''(/,5) in terms of fermion and boson oper- 
ators /, 6, respectively. For instance, in the very pop- 
ular slave-boson representation one has Xf"' — fi^ab\ 
with the constraint b\hi -I- Ylia=i fiafi,<^ ~ 1 where fa 
destroys a spinon with spin a and w creates a holon. 
Since fa- , b obey the canonical Fermi and Bose commu- 
tation rules it is, at first glance, very convenient to use 
this approach since for these operators well-elaborated 
many-body techniques (Feynman diagrams) exist. How- 
ever, the constraint makes the spinon and the slave boson 
strongly coupled, and a naive decoupling of the equations 
of motion is risky, especially in dimensions d ^ 2, some- 
times giving fictitious results. Therefore, the present au- 
thor and his collaborators decided to work directly with 
the Hubbard operators, thus always keeping the compos- 
ite object (bound spinon and holon in the language of 
SB). However, the central question in this approach is - 
how to study the quasiparticle dynamic in absence of the 
canonical Fermi and Bose commutation rules? The way 
out was found in the very elegant Baym-Kadanoff tech- 
nique in obtaining Dyson equations, the method which 
is independent on the operator algebra - see more in [7| 
and references therein. Applied to strong correlations we 
call it the X-method. It turns out that the method is very 
efficient and one can formulate very elegant and control- 
lable 1/N expansion which also allows us to study the 
EPI interaction in a consistent way. We stress that until 
now there has been no transparent and efficient way to 
treat the EPI within the slave-boson method. 

The basic idea behind the X-method is that the 
Dyson's equation for the electron Green's function can be 
effectively obtained by introducing the external potential 
M°'i°'^(l) (source) with the Hamiltonian Hs 



H,dT 



E 



M'"i'"^(l)X'"i'"=(l)dl 



^(l)X'^i'"^(l) 



(13) 



Jy (Si • Sj 



-fliflj ) + H3. 



(12) 



The first term describes electron hopping by ensur- 
ing that the double occupancy is excluded. The sec- 
ond term describes the Heisenberg-like exchange energy 



Here, 1 = (1,t) and J{..)dl = J{- )dTj2i and r is the 
Matsubara time. In the following, integration over the 
bared variables (1,2..) and a summation over bared spin 
variable (a..) is understood. The sources ^'^'^'^^(l) are 
useful in generating higher correlation functions entering 
the self-energy. The electronic Green's function is defined 
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by 0, M 

-{f fs'i'0'"i(l)X^^0(2))} 

Qa^a2 3) = !^ ^, (14) 

(TS) 

where T is the time-ordering operator and S = 
Texp{— J Hs{l) dl}. We define the quasiparticle Green's 
g^j-^j ^Yie Hubbard spectral weight Q'^^'^^ 



— 1,(T2'72 



(2) 



g'"i'"^(l) = (5'"i'"^(X°"(l)} + (X'^^'^^il)). (15) 
It is useful to define the vertex functions 

J^g-l,<Tl<T2(l^2) 



(16) 



which in what follows play very important role in the 
quasiparticle dynamics. It turns out that the self-energy 
(note = g^^ — Eg - 0) can be expressed via 
the charge vertex 7^,(1, 2; 3) = 7^^(1,2; 3) and the spin 
vertex 7^(1, 2; 3) = 7^^(1,2; 3). In the paramagnetic 



state one has .g'^^'^^ 



Tia2 = S^'^'Q and 



5^.(l-2) = -^i^Q(l) 
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Ml ~ 1) 



N 



^o(l-l) 
iV 



5(1 -2)7,(2, 2; 1) 



+^2(1, 1,3) 

TV 



3(1-2)7,(2,2;3) + S]q(1-2) 



(17) 



where t2{l, 2, 3) = 6{1 - 2)io(l - 3) - (5(1 - 3) Jo(l - 2). 
The notation to{l — 2) (and Jo{l — 2)) means to{l — 2) — 
^o.iij2^{'''i ~ ''"2)- The first two terms in Eq. (fTTjl repre- 
sent the effective (reduced) kinetic energy of quasiparti- 
cles in the lower Hubbard band. As we shall see below 
they give rise to the quasiparticle band narrowing and to 
the shift of the band center, respectively. The third and 
fourth terms describe the kinematic and dynamic interac- 
tion of quasiparticles with charge and spin fluctuations, 
respectively, while the term I]q(1,2) takes into account 
the counterflow of surrounding quasiparticles which takes 
place because of the local constraint (absence of double 
occupancy). E depends on the vertex functions 7c(l, 2; 3) 
and 7j;(l, 2; 3) and it does not contain a small expansion 
parameter, like the interaction energy in weakly interact- 
ing systems, because the hopping parameter t describes 
at the same time the kinetic energy and kinematic in- 
teraction of quasiparticles. This means that there is no 



controllable perturbation technique for Eg, due to the 
lack of a small parameter. In the past various decou- 
pling procedures and mean-field like techniques were uti- 
lized, such as the steepest descent method in the path 
integral technique, etc. However, to the author's best 
knowledge, all decoupling techniques are not only non- 
systematic and uncontrollable, but moreover they suffer 
from inability to extract the coherent part of the quasi- 
particle spectrum. As we shall argue in what follows the 
coherent quasiparticle band-width (in the case J = 0) is 
proportional to 5, i.e. Wqp ^ 6 ■ t. All decoupling tech- 
nique obtain the band- width Wd ^ {1 — a ■ 5) • t [53l| . 
with a ■ S finite in the limit 6—^0. So obtained Wd is 
in fact the total band-width for coherent and incoherent 
quasiparticle motions, i.e. Wd = Wqp + Wine- Thus the 
decoupling techniques are inappropriate in studying the 
coherent quasiparticle motion. 

What is then the advantage of the X-method and the 
self-energy in Eq. IjlTII ? This method allows formulation 
of a controllable 1 /N expansion for Eg by including also 
the EPI, as shown below. For that purpose it is necessary 
to generalize the local constraint condition [s^l 
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(18) 



where N/2 replaces the unity in Eq. Hll|l . Only in this 
case the 1/N expansion is systematic. It is apparent from 
Eq. I^HJ that for = 2 it coincides with Eq. and 
has the meaning that at most half of all spin states at a 
given lattice site can be occupied. The particle spectral 
function A{\i,Lu) = —ImG{\<i^uj) /it must obey the gen- 
eralized Hubbard sum rule which respects the new local 
constraint in Eq. H18|l 



dLoA{]<i, uj) 



1 + {N - 1)5 



The N > 2 generalization of the local constraint allows 
us to make a controllable 1/N expansion of the self- 
energy with respect to the small quantity 1/A^ (when 
A^ ^ 1). Physically this procedure means that we are se- 
lecting a class of diagrams in the self-energy and response 
functions which might be important in some parameter 
regime. In [33|, [Ij, it was shown that there is a 
systematic 1/A^ expansion for vertices 7, and 7^ as well 
as for other functions 



9 = 90 +91/N + 

Q = Nqa + Qi + Q2/N.. 

E = Eo + Ei/Ar + ... 



(19) 



etc. 



1. The main results of the X-method 
ms in lea ding 



The first nontrivial terms in lea ding O(l)-order are Eq 
and go - for details see in , [s^l , [Ul , [Sj] , which in fact 
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describe the coherent part of the quasiparticle self-energy 
and the Green's function 

go{k,uj) = Go{k,uj)/Qo = , 

u - (eo(k) - n 

eo(k) = ec - goio(k) - — ^ Jo(k + p)nF(p), 

^ p 

^ p 

Qo = < X^^ >= Nqo = N^. (20) 

where iV^ is the number of lattice site, eo(k) is the quasi- 
particle energy for the coherent motion and is the 
level shift. Here, to(k) and Jo(k) are Fourier trans- 
forms of to^ij and Jo.ij, respectively. For t — t' one 
has io(k) — 2to{coskx + cos ky) + 2t' cos cos ky and 
Jo(k) ~ 2Jo{coska: + COS ky), since the lattice constant is 
let a = 1. Since Qo{^ Nd) this means that the quasipar- 
ticle residuum vanishes at vanishing of doping {S —^ 0). 
This is physically plausible since for (5 ^ the coherent 
motion in the t-J model is blocked. The chemical poten- 
tial /i is obtained from the condition 1 — 5 = 2 "f(p) 
which gives large Fermi surface, see discussion below. 
Note that Eg in Eq. (HTjl is of the 0{1/N) order. 

We summarize the main results obtained by the X- 
method to the leading 0(1) order and compare them 
with the corresponding ones in the slave boson (SB)- 
method [^2: (1) up to the 0(1) order, 5o(k, w) 

describes the coherent motion of quasiparticles in the 
metallic state, whose contribution to the total spec- 
tral weight of the particle Green's function Gq (k, uj) is 
Qo = NS/2, i.e. G'o(k, w) = (5offo(k, t^^) to the leading 
order, so we dealing with a Landau-Fermi liquid. This 
means that the quasiparticle residuum (and the metallic 
state) vanishes for zero doping ^ = and the system is 
in the Mott insulating state. The quasiparticle energy is 
dominated by the exchange parameter if Jq > 6to, i.e. 
for very low doping, since in cuprates (Jo/^o) < 1/3. For 
Jo = there is band narrowing since 6 < 1 and the quasi- 
particle band- width is proportional to S, i.e. W = z-S-tQ. 
These results are identical to the corresponding ones in 
the SB method. This does not mean that the next lead- 
ing terms coincide in both methods. This is actually 
not the case for a number of quantities (response func- 
tions). This is understandable since the X- method keeps 
the composite object (correlated motion of spinon and 
holon), while in the SB method they are decoupled. Only 
after invoking the existence of gauge fields (also fictitious 
particles) one can keep spinon and holon together in the 
SB model, if at all. (2) The X-method respects the lo- 
cal constraint at each lattice site in every step of calcu- 
lations, which guarantees the correct study of response 
functions. (3) In the (important) pa p er |5 4l| - which is 
based on the theory elaborated in [s^l, |3lj. |33|. it was 
shown that there is large discrepancy in approximate cal- 
culations within the SB and X methods. For instance, 
in the superconducting state the anomalous self-energy 



(which is of 0(l/A^)-order in the 1/N expansion) of the 
X- and SB methods differ substantially. As a result, the 
SB method erroneously predicts superconductivity due 
to the kinematic interaction in the t-J model (for J = 0) 
with large T^, while the X- method gives extremely small 
Tc (« 0) The reason for this discrepancy is that 

the calculations done within the SB method miss a class 
of compensating diagrams, which are in contrast taken 
into account automatically within the X method. So, 
although the two approaches yield some similar results 
in leading 0(1) order their implementation in the next 
leading 0(1/A^) order show that they are indeed different 
- see discussion below. Note that the 1/N expansion in 
the X-method is well-defined and transparent. (4) The 
renormalization of the EPI coupling constant by strong 
correlations is different in the two approaches even in 
the large A^-limit, as shown below; (5) The optical con- 
ductivity cri(ti;,q = 0) = ai{uj) and the optics sum-rule 
exhibit very interesting behavior as a function of doping 
concentration 

J duja{uj) = -d ■ Ne^a^ec- (21) 

It is seen that the optics sum-rule is proportional to the 
number of holes 5, instead of n = 1 — i5 as it would be in 
the weakly interacting case . Here a is the lattice con- 
stant. (6) Note, that the volume below the Fermi surface 
in the case of strong correlations scales with n = 1 — 5, i.e. 
the Fermi surface is large like in the conventional Fermi 
liquid. The above analysis clearly shows significant dif- 
ference ((5) vs (6)) in response functions of strongly cor- 
related systems and the canonical Landau-Fermi liquid. 



IV. THE FORWARD SCATTERING PEAK IN 
THE EPI OF CUPRATES 

The EPI coupling constant in LTSC materials is cal- 
culated by the local-density functional (LDA) method. 
The latter is suitable for ground state properties of crys- 
tals and it is based on an effective electronic crystal po- 
tential Vg . Since in principle Vg differs significantly from 
So(k, UJ = 0) (the many body effective potential - see 0) 
then the LDA coupling constant g^j^p^^ can be also very 
different from the real coupling constant gEP- Strictly 
speaking the EPI does not have meaning in the LDA 
method - see discussion in 0, because the latter treats 
ground state properties of materials, while the EPI is due 
to excited states and inelastic processes in the system. 
We shall not deal with this problem here - see 0- 



A. LDA method for the EPI in cuprates 

The LDA method considers electrons in the ground 
state (there is a generalization to finite T), whose energy 
can be calculated by knowing the spectrum {efc} of the 
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Kohn-Sham (Schrodinger like) equation 



[;^ + K,(r)]^,(r)=efcV,(r) 
2m 



(22) 



which depends on the effective one-particle potential 

Vg{v) = K^(r) + Vh{v) + Vxc{r). (23) 

Here, Vei is the electron-lattice potential, Vh is the 
Hartree term and Vxc describes exchange-correlation ef- 
fects - see 0- Because the EPI depends on the excited 
states (above the ground state) of the system, this means 
that in principle the LDA method can not describe it - 
see 0- However, by using an analogy with the micro- 
scopic Migdal-Eliashberg theory one can define the EPI 
coupling constant g^^^s) = gVc/s also in the LDA theory, 
see H- It reads 
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where n means summation over the lattice sites, a — 
x,y,z and the wave function V'ki the solution of the 
Kohn-Sham equation. Formally one has dVg/SlEin — 
^LDA^e^^Vei- Evcn in such a simplified approach it is 
difhcuh to calculate gl^if,^^ = fff.^^^^ + ffS^'"'' because 
it contains the short-range (local) coupling 

{Ylrn I I Yi.m') (25) 
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fffr-'(k,k') 



with AZ = 1, and the long-range coupling 

g2T'{k,k')^{Yi„,\{K-IllX\Yi,. 



(26) 



with Al ~ 0. In most calculations the local term g^n^"^ 
is calculated only, which is justified in simple metals only 
but not in the HTSC oxides. In cuprates the local term 
gives a very small EPI coupling x^'^^"^^ ~ 0.1, which 
is apparently much smaller than the experimental value 
A > 1, and which gives rise to a pessimistically small 
Tc 113. The small X^'^'^^ was also one of the reasons 
for discarding the EPI as pairing mechanism in HTSC 
oxides. At the beginning of the HTSC era the electron- 
phonon spectral function a^F(a;) was calculated for the 
case La2-xSrxCuOi in 's^l by using the first-principles 
band structure calculations and the nonorthogonal tight- 
binding theory of lattice dynamics. The value A = 2.6 
was obtained and with the assumed /x* — 0.13 this gave 
Tc = 36 K. However, these calculations also predicted 
a lattice instability for the oxygen breathing mode in 
Lai.^^SrQ i^CuO/^ that has never been observed. More- 
over, the same method was applied to Y Ba2Cu^Oj in 
[ssj l where it was found A = 0.5, which at best provides 
Tc = (19 - 30) K. In fact the calculations in ^ did 



not take into account the Madelung coupling (i.e. they 
neglected the matrix elements with Al — Q). 

However, because of the weak screening of the ionic 
(long-range) Madelung coupling in HTSC oxides - es- 
pecially for vibrations along the c-axis - it is necessary 
to include the nonlocal term This task was 

achieved within the LDA approach by the Pickett's group 
psj . where the EPI coupling for La2~xMxCuOi is cal- 
culated within the frozen phonon method for q 
modes. By extrapolating the result for q = to fi- 
nite q they obtained A = 1.37 and wiog ~ 400 K, and 
for ^* = 0.1 they predicted that = K (T^ « 
cjiog exp{-l/[(A/(l -I- A)) - n*]})- The justification of 
this procedure was questioned in [? ] were smaller A 
was obtained. For more details see Ref. Q and ref- 
erences therein. Next, some calculations, based on the 
tight-binding parametrization of the band structure in 
Y Ba2Cu?, Oj, g ave rather large EPI coupling A « 2 and 
Tc = 90 if 

We point out, that model calculations which take into 
account the l ong -range ionic Madelung potential appro- 
priately |5^ , |57j| also gave rather large coupling constant 
A ^ 2, which additionally hints to the importance of long- 
range forces in the EPI - see and references therein.. 

Since in HTSC oxides the plasma frequency along the 
c-axis, Wpj, is of the order (or even less) of some charac- 
teristic c-axis vibration modes, it is necessary to include 
the nonadiabatic effects in the EPI coupling constant, 
i.e. its frequency dependence ga.n ^ This 
non-adiabaticity was partly accounted for by the Falter 
group ^% by calculating the electronic dielectric func- 
tion along the c-axis ecc(k, uj) in the RPA approximation. 
The result was that ga,n increased appreciably beyond its 
(well screened) metallic part, which gave a large increase 
of the EPI coupling not only in the phonon modes but 
also in the plasmon one. This question deserves much 
more attention than what it received in the past. 

The electron band structure and the Fermi surface can 
be satisfactory described by the tight binding model, 
which is characterized with various hopping parameters 
tij and the local atomic (ionic) levels e;. Accordingly, 
there are there two kinds of the EPI coupling. The co- 
valent part of the EPI is due to strong covalency of the 
Cu and O orbital in the CuO^ planes. In this case, the 
EPI coupling constant is characterized by the parame- 
ter ("field") E^'"' ~ dtp_d/dR ~ qcovtp-d, where 
is the hopping integral between Cu^dx^-y'^) and 0{px.y) 
orbital and the length q~^^ characterizes the spacial ex- 
ponential fall-off of the hopping integral tp^d- As the 
phonon Raman scattering shows, the covalent EPI is un- 
able to explain the strong phonon renormalization (the 
self-energy features) in the Big mode in YBa2Cu30j by 
superconductivity, since in this mode the 0-ions vibrate 
along the c — axis in opposite directions and for this 
mode dtp^d/dR is zero in the first order in the phonon 
displacement. Therefore, the EPI in this mode is cer- 
tainly due to the ionic contribution which comes from 
the change in Madelung energy, as it was first proposed 
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in [s^, • Namely, the Madelung interaction creates 
an electric field perpendicular to the CUO2 planes, which 
is due to the surrounding ions that form an asymmetric 
environment. In that case the site energies e° contain the 
matrix element e*°" — \ V{y) \ 1/;^), where | -0^) is the 
atomic wave function at the i-th site, while the poten- 
tial V{r) steams from the surrounding ions. In the sim- 
ple and transition metals the surrounding ions are well 
screened and therefore the change of e'°" in the presence 
of phonons is negligible, in contrast to cuprates which are 
almost ionic compounds, in particular along the c-axis 
where the change of e'°" is appreciable and characterized 
by the field strength — Vivi/dn- Here, Vm is the 

characteristic Madelung potential due to the surround- 
ing ions and d„ is the distance between the neighboring 
ions. Immediately after the discovery of cuprates it was 
assumed in many papers [58j that the covalent part dom- 
inates the EPI in these materials. However, the calcula- 
tion that considered only the covalent effects 58] gave a 
rather small (-10-20 K in YBCO, and 20-30 K in 
Lai,^^SrQ i^CuO/i). It turns out that in cuprates the op- 
posite inequality E"^"^ ^ E'^°'" is realized for most c-axis 
phonon modes in spite of the fact that qcov > ^/dn - 
see more in and references therein. This is supported 
by detailed theoretical studies of YBCO |H3| , where the 
change in the ionic Madelung energy due to the out of 
plane oxygen vibration in the Big mode is calculated . 
Similar to YBCO, the large superconductivity-induced 
phonon self-energy effects in HgBa2Ca^Cu/iOioj^x and 
in {Cu,C)Ba2Ca3Cu40io-f-x for the Aig modes are also 
due to the ionic (Madelung) coupling. In these modes 
oxygen ions move also along the c-axis and the ionicity 
of the structure is involved in the EPI. This type of long- 
range EPI is absent in usual isotropic metals [LTSC 
superconductors), where the large Coulomb screening 
makes EPI local. Similar ideas were recently incorpo- 
rated into the Eliashberg equations in |5^. The weak 
screening along the c-axis, which is due to the very small 
hopping integral for carrier motion, is reflected in the 
very small plasma frequency uip along this axis. Since 
for some optical phonon modes one has Uph > lu^"^ , nona- 
diabatic effects in the screening are important. The latter 
can give rise to much larger EPI coupling constant for this 
modes |23| . These ideas are also supported by the recent 
ARPES measurements on the single-layer Bi2Sr2CuO(, 
(Bi2201) where the coupling of electrons with phonons in 
the region 30—60 meV in the overdoped compound is sig- 
nificantly reduced compared to the optimally doped case. 
This can be explained by the larger electronic screening 
along the c-axis in overdoped compound. 

To summarize, ARPES, electron and phonon Raman 
scattering measurements in the normal and supercon- 
ducting state of cuprates gave the following important 
results: (a) phonons interact strongly with the electronic 
continuum, i.e. EPI is substantial; {b) the ionic contri- 
bution (the Madelung energy) to EPI interaction for c- 
axis phonon modes gives substantial contribution to the 
(large) EPI coupling constant (A > 1). 



B. Renormalization of the electron-phonon 
interaction by strong correlations 

Based on the above discussion the minimal (toy) model 
Hamiltonian for HTSC cuprates contains besides the t-J 
terms also the electron-phonon interaction 

H = HtJ + ^ (^^a,i^i" + Hph 
i.<7 

where Hepi H'^p are the ionic and covalent contribu- 
tion to the EPI, respectively. We consider first the ionic 
term where ^^(ul,^) describes the change of the atomic 
energy e°_j due to the change of long-range Madelung 
energy in the presence of phonon displacements vll^k, 

of 

other atoms. L and k enumerate unit lattice vectors and 
atoms in the unit cell, respectively. Note that the theory 
is formulated for the general nonlinear form of $i(ui „) 
dependence, and the following analysis holds in principle 
also for an anharmonic EPI. The term proportional to 
(Xf ) in Eq. (|77|) is introduced in order to have = 
in the equilibrium state. We stress that the X-method 
allows controllable calculation of the covalent contribu- 
tion In contrast, the treatment of EPI by the SB 
method is complicated and not well defined, giving some- 
times wrong results - see more in 7J . 

After technically lengthy calculations, which are per- 
formed in [33, |3l|. the expression for the ionic contribu- 
tion to the EPI part (frequency-dependent) of the self- 
energy reads 

t!^;^T\i,2) = -FBp(i, 2)7,(1, 3; 1)50(3, 4)7,(4, 2; 2) 
Vep{1,2) = e-\l ~ l)V^p{l - 2)e-\2 - 2) 



yi!p(l-2) = -(T|.(l)$(2) (28) 

The propagator Vgp(l — 2) (which includes also the cou- 
pling constant) of the bare EPI comprises in principle 
also the anharmonic contribution. From Eq. (|28|) it is 
seen that in strongly correlated systems the ionic part of 
the EPI is proportional to the square of the three-point 
charge vertex 7,(1, 2; 3) (due to correlations). The self- 
energy is given by 

/"OO 

4^]^") (k, uj) = / dn{a'F{\^, k', 17))k'i?(c^, ^) 
Jo 
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where the (momentum-dependent) Ehashberg spectral 
function a^F(k, k', uj) is defined below, nsi^) is the Bose 
function and ijj{x) is the di-gamma function - 0, |30j |. 



C. Forward scattering peak in the cliarge vertex 7^ 

The three-point charge vertex 7^(1, 2; 3) plays an im- 
portant role in renormalization of all charge processes, 
such as EPI in Eq. (|28|l . Coulomb scattering and the 
scattering on non-magnetic impurities. 7^(1, 2; 3) (and 
its Fourier transform 7^(k, g), q = (q, ign), = 2'!mT,) 
was calculated in the t — J model in leading 0(1) order 
|30| , [U , 112 , ■ It is the solution of the following linear 
integral equation 

7,(l,2;3) = <5(l-2)(5(l-3) 



+to{l - 2)50(1 - l).go(2 - l+)7c(l, 2; 3) 



+S{1 - 2)to(l - 1)50(1 - 2)50(3 - l)7c(2, 3; 3) 



Jo(l - 2)50(1 -l))5o(2- 2)7,(1, 2; 3) 



(30) 



The (charge) vertex function 7c(k, q) describes a specific 
screening of the charge potential due to strong correla- 
tions. In the presence of an external perturbation (u) 
there is a change of the band-width, as well as the change 
of the local chemical potential, which comes from the sup- 
pression of double occupancy. These processes are con- 
tained in the second and the third term of Eq. 1301) . The 
central result is that for momenta k laying at (and near) 
the Fermi surface, 7,o('^i*li^ — 0) ^ strong forward 
scattering peak at q = 0, which is very pronounced for 
lower doping 5 (^ 1). On the other hand, the backward 
(at large q) scattering is substantially suppressed, see 
Fig. Such a behavior of the vertex function means 
that a quasiparticle moving in the strongly correlated 
medium digs up a giant correlation hole with the radius 
^ch ~ '^/^' where a is the lattice constant. The effect of 
the J-term is small, since 7^, is determined by charge fluc- 
tuations which are spread over the whole bare bandwidth 
Wb zto » Jo- 

We stress that 7,9 Fig- calculated in the adia- 
batic approximation, i.e. for w = in 7,0(^1 Qi '-^)- How- 
ever, in the nonadiabatic regime cj > q • v^(q) the vertex 
function reaches its maximal value 7,(ki?, q = 0,0;) = 1, 
see Fig. 0j3. This nonadiabaticity occurs whenever the 
phase velocity (w / q) of charge excitations is larger than 
the Fermi velocity. This effect is realized in case of high 



OA 



0.3 



cr 

0.2 

X 

CD 
> 

0.0 



-0.1 



-0.2 





J=0.2, t'=0. 




q=(q-q) 


\ \ 


5:^:0.50 






I \ 5=0 20 




^ \- 








\ 




\ 5=0.10 









0.0 



1.0 2.0 

momentum aq 



3.0 



3 

cr 0.5 

LL 

^ 0.4 













(b) 















aq 



aq 



FIG. 4: (a) Adiabatic {u> = 0) vertex function 7(kF, q) of the 
t-J model as a function of the momentum aq with q = {q, q) 
for three different doping levels S. From Ref. fs^. (b) Non- 
adiabatic {u 7^ 0) vertex function ■y{\i.F ,ci,uj) schematically 
with q^ — lu/vf- 



energy phonons. This means that due to nonadiabatic 
effects strong correlations do not suppress the EPI cou- 
pling with optical phonons at g = 0, in contrast to a 
number of claims |63| that they do. 

Since the real physics is characterized by TV = 2 one 
can ask - how reliable are the calculations of 7,(ki?,q) 
by the 1/A^ expansion method and in the limit N 00. 
Here let us mention that VL (Ginzburg) used to para- 
phrase a Landau's remark "that the many-body theory 
is an experimental science" . This means that the cor- 
rectness of the theory is tested by experiments. In this 
respect, the exact diagonalization of charge correlation 
function A^(q, uj) in the t-J model shows clearly that 
low-energy charge scattering processes at large momenta 
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I q jw 2k p are strongly suppressed compared to scatter- 
ing at small transferred momenta (| q |^ 2kp). This 
unambiguo usly confirms the results obtained by the X- 
method in '30] , 'sT], "32] on the suppression of backward 
scattering in the vertex. Second, very recent Monte Carlo 
(numerical) calculations of the vertex function "/^(kp, q) 
in the Hubbard model with finite U 62] show clear de- 
velopment of the forward scattering peak in 7^(ki?,q), 
that is more pronounced for larger U. This result defi- 
nitely confirms the predictions of appearance of the for- 
ward scattering peak in ^^(kp,q) in Rcfs. [s^l, [13, [s^ . 
• For analytical theories numerical calculations in some 
models play the role of an experiment, then the above 
theory of strong correlations and EPI is in that sense 
confirmed experimentally in accord with the Landau re- 
mark. 

Finally, it is possible to calculate the contribution 
of covalent EPI coupling H'^p to the self-energy. In 
this case, the vertex renormalization by strong corre- 
lations is determined by the four-point vertex function 
7((1, 2; 3, 4). The latter is an important ingredient of EPI 
not only in HTSC cuprates but also in heavy fermion 
systems. The calculations in show that 7^ differs 
from that obtained by the mean-field approximation SB 
model, what implies that the previous studies of EPI in 
heavy fermion materials, in which the "covalent" cou- 
pling ~ 7( dominates, should be reconsidered by using 
the correct form of 74. 



D. Pairing and transport EPI coupling constants 

Depending on the symmetry of superconducting order 
parameter A(k, w) (s— , d-wave pairing) various averages 
(over the Fermi surface) of a^i^(k, k', uj) enter the Eliash- 
berg equations. Assuming that the superconducting or- 
der parameter transforms according to the representa- 
tion Ti (i = 1,3,5) of the point group of a square 
lattice (in the Cu02 planes), the appropriate symmetry- 
projected spectral function is given by 

a2i^.(k,k ,tu) = I go,«c.(k,k-r,k',i^) p X 

X 5{u - o.,(k-r,k')) I 7c(k, k-T,k') |2 A(j), (31) 

where k and k' are the momenta on the Fermi line 
in the irreducible Brillouin zone which is equal to 1/8 
of the total Brillouin zone. Here, (7o,scr(k, p, i^) is the 
EPI coupling constant for the z^-the mode, where the 
screening by long-range Coulomb interaction is included, 
i.e. 5o,scr(k, p, i^) = 5o(k,p, i^)/ee(p). The density 
of states 7Vsc(0) is renormalized by strong correlations, 
where N^cifi) = No{0)/qo and go = S/2 in the t-t' 
model {J — Q). In the t — J model Nsc{0) has an- 
other form which does not diverge for ^ — > 0, but one 
has iVsc(0)(~ 1/Jo) > -^o(O), where the bare density 



of states iVo(O) is calculated, for instance by the LDA 
scheme. Tj , where j = 1, ..8, denotes the eight point- 
group transformations that form the symmetry group of 
the square lattice. This group has five irreducible repre- 
sentations which we distinguish by the label i — 1,2, ...5. 
In the following the representations i — 1 and i — 3 
, which correspond to the s— and d-wave symmetry of 
the full rotation group, respectively, will be of impor- 
tance. Di{j) is the representation matrix of the j-th 
transformation for the representation i. For each sym- 
metry one obtains the corresponding spectral function 
a'^Fi{uj) = {{a'^Fi{k,k ,oj))^)^,, which (in the first ap- 
proximation) determines the transition temperature for 
the order parameter with the symmetry F^. In the case 
i — 3 the electron-phonon spectral function 0^-^3(0;) in 
the d- channel is responsible for d-wave superconductiv- 
ity represented by the irreducible representation Fa (or 
sometimes labelled as Big). 

Performing similar calculations (as above) for the 
phonon-limited resistivity one finds that the latter is 
related to the transport spectral function a^Ftrii^) = 
((a2F(k,k',a;)[v(k) - v(k')]2)k)k'/2((v2(k))k)k', where 
v(k) is the Fermi velocity. The effect of strong correla- 
tions on EPI was discussed in 13^ and more extensively 
in 31], 32] within the model where the phonon frequen- 
cies a;(k — k ) and (7o,scr(k, p. A) are weakly momentum 
dependent. In order to illustrate the effect of strong cor- 
relations on ct^Fi{uj) we consider the latter functions at 
zero frequency (w = 0) which are then reduced to the (so 
called) "enhancement" functions 

= ^^E«l 7e(k,k-T,k') Wh,D..{3). (32) 

Note that in the case J = one has Nsc{0)/No{0) = q^^, 
where qq is related to the doping concentration, i.e. 
qq = (5/2. Similarly, the correlation effects on the re- 
sistivity /3(T)(~ Atr) renormalize the transport coupling 
constant At^ (similarly to A^ it is defined via a^Ftritt')). 
Note that for quasiparticles with an isotropic band the 
absence of correlations implies that Ai = Atr = 1, A; = 
for i > 1. The averages in A 1, A3 and At^, shown in Fig.|Sl 
were performed numerically in [sij by using a realistic 
anisotropic band dispersion in the t-t'-J model and the 
corresponding charge vertex. The three curves are multi- 
plied by a common factor so that Ai approaches 1 in the 
empty-band limit 5 — > 1, when strong correlations are 
absent. Note that Tc in the weak-coupling limit and in 
the i-th channel scales like Ti*-* « (lu) exp(— l/(Ai — ^*), 
where /z* is the Coulomb pseudopotential in the i-th 
channel and (lo) is the averaged phonon frequency. 

Several interesting results, which are shown in Fig. |S1 
should be stressed. First, in the empty-band limit S ~> 1 
the d-wave coupling constant A3 is much smaller than 
the s-wave coupling constant Ai, i.e. A3 ^ Ai. Further- 
more, the totally symmetric function Ai decreases with 
decreased doping. Second, in both models Ai and A3 ap- 
proach one another at some small doping (5 w 0.1 — 0.2, 
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where Ai « A3 but still Ai > A3. By taking into account 
residual Coulomb repulsion of quasiparticles which usu- 
ally have /i^ «C /^^ one finds that the s-wave supercon- 
ductivity (which is governed by the coupling constant 
Ai) is strongly suppressed by /i*, while the d-wave su- 
perconductivity (governed by A3) is only weakly affected 
by fi^- In that case the d-wave superconductivity due 
to EPI becomes more stable than the s-wave supercon- 
ductivity at sufficiently small doping S, i.e. Tc"^^ > Tc'^K 
Interference experiments |6l| show that this occurs in 
underdoped, optimally doped and overdoped cupratcs. 
This means, that EPI is responsible for the strength of 
pairing in cuprates, while the d-wave superconductiv- 
ity is triggered by the residual Coulomb interaction (in- 
cluding also spin- fluctuation scattering). Third, in the 
nonadiabatic regime when the phonon frequency fulfills 
the condition ojph > q • (q) the enhancement func- 
tion 7^(kp, p,ix'p;i)/(3'o is substantially larger compared 
to the adiabatic one, which is shown in Fig. |S1 This 
means that the strength of EPI coupling is differently af- 
fected by strong correlations for different phonons. For a 
given frequency the coupling to phonons with momenta 
p < Pc — i^/vp enhanced, while the coupling to those 
with p > pc = Lo/vp is substantially reduced due to sup- 
pression of backward scattering by strong correlations. 
Fourth, the transport coupling constant Atr (not prop- 
erly normalized in Fig. see the correction in |32l|) is 
reduced in the presence of strong correlations, especially 
for lower doping where one has approximately Atr < A/3 
|32|. This is an important result since it resolves the 
experimental puzzle - that Atr (which enters resistivity 
p{T) ~ XtrT) is much smaller than the coupling constant 
A (which enters the self-energy S and Tc), i.e. Atr << A. 
To remind the reader, Atr ~ A is realized in almost all 
LTSC. One of the important differences between LTSC 
and HTSC in cuprates lies in strong correlations present 
in cuprates which causes the forward scattering peak in 
charge scattering processes. Fifth, the forward scattering 
peak in EPI of strongly correlated systems is a general 
phenomenon affecting electronic coupling to all phonons. 
This means that the bare EPI coupling constant for each 
phono-mode (v) must be multiplied by the vertex func- 
tion, i.e. 5°(k,q) ^ 7(k, q)5^(k, q). 

As we already said, Monte Carlo (numerical) calcula- 
tions of the Hubbard model at finite U - performed by 
Scalapino's group [63|, show that the forward scattering 
peak in the EPI coupling constant (and the charge ver- 
tex) develops as U is increased. These numerical results 
prove the essential correctness of EPI theory based on 
the X method. The latter effect is more pronounced at 
lower doping. Results similar to Monte Carlo ones were 
obtained recently within the framework of Riickenstein- 
Kotliar (four slave-boson) model [63l |. 

We stress that, contrary to the X-method where the 
systematic 1/N calculations of the EPI self-energy is 
uniquely done, this is still a problem for the SB method 
where an 1/iV expansion of the partition function Z(T, fi) 
is usually performed |64j . The expression for the vertex 
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function in the SB method is different from that in the 
X method 30], 31]. It seems that the existing SB treat- 
ment of EPI omits a class of diagrams which causes an 
incorrect dependence of the coupling constant A on dop- 
ing. As a result, the vertex function in the SB approach 
is peaked not at g = but at some finite (/max, where 
gmax = only for doping 5 = 0. 
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V. NOVEL EFFECTS DUE TO THE FORWARD 
SCATTERING PEAK IN EPI 

There are a number of effects which are predicted by 
the theory with FSP in EPI and in nonmagnetic impurity 
scattering - the FSP theory. We have already explained 
the effects of the forward scattering peak on the EPI. We 
discuss briefly some other predictions of the FSP theory 
including aspects not covered in 0- 



A. Theory of ARPES self-energy in cuprates 

Recent ARPES experiments of the Shen group [20| . 
H^i inSli [H^ gave additional evidence that phonons are 
involved in pairing mechanism of cuprates. Furthermore, 
the clear isotope effect was found in the real part of the 
ARPES self-energy Rel](k, w) "2^1, which confirms the 
importance of EPI interaction in cuprates. 



assume that the former case is realized, although this 
is not crucial because ARPES indicates that electron- 
phonon coupling is much larger than the Coulomb cou- 
pling, i.e. Xpfi » Ac; (iv) The scattering potential from 
non-magnetic impurities has pronounced forward scat- 
terin g pe ak, which is also due to strong correlations |3n| | , 
[3ll |. |32||. The latter is characterized by two rates 7i(2)- 
The case = 72 mimics the extreme forward impurity 
scattering. In this case, d-wave pairing is unaffected by 
impurities - the Anderson theorem for unconventional 
pairing. On the other hand the case 72 = describes the 
isotropic exchange scattering - see discussion in • 
The Green's function is given by 



Gk 



(33) 



where k = (k, lo) . In the FSP theory the equations for 
Cbk and Afc are local on the Fermi surface, i.e. in k space 
and Wi 



1. ARPES kvnk and non-shift puzzle 

Recently, it was reported that in all metallic cuprates 
there is a kink in the quasiparticle dispersion (^(^k) the 
nodal direction (along the (0,0) — {tt,7t) line) at around 



,(70) 



Sec.II.F 



(60 



70) meV |20|, see the property (Gn) in 
and Fig. |21 In the superconducting state this 
kink is not shifted (the property (3s) in Sec. II. F) con- 
trary to the standard Eliashberg theory - we call this 
the ARPES non-shift puzzle. In the case of the anti- 
nodal point (7r,0) there is a kink at 40 meV in the nor- 
mal state which is shifted by the maximal gap in the 
superconducting state (the property (4^) 23], see Fig.|31 
This means that any theory which aspires to explain the 
pairing in cuprates must solve the non-shift puzzle. It 
is worth mentioning that in spite of the apparent exis- 
tence of the non-shift puzzle since 2001, practically there 
have been no publications that treat this problem ex- 
cept the paper l37| where a plausible theory was given. 
The approach in |33 is based on the FSP theory for EPI 
and other charge scattering processes, which is discussed 
above. In order to explain the non-shift puzzle by the 
FSP theory the following simplifications are made: (i) 
electron-phonon interaction is dominant in HTSC, and 
its spectral function Q!^F(k, k', fi) « a'^ F {if , ip' , Q) {ip 
is the angle on the Fermi surface) has a pronounced 
forward scattering peak due to strong correlations. Its 
width is very narrow | k — k' kp even in over- 

doped systems [s^, Hj, [l^. To the leading order, 
one can put a'^F{(p,(p\fl) ^ 5{ip — ip'); (ii) for sim- 
plicity the slightly broadened Einstein spectrum is as- 
sumed; (Hi) the dynamic part (beyond the Hartree-Fock) 
of the Coulomb interaction is characterized by the spec- 
tral function S'c(k,k',0). The ARPES non-shift puzzle 
implies that Sc is either peaked at small transfer mo- 
menta I k — k' I, or it is so small that the shift is weak 
and beyond the experimental resolution of ARPES. We 
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where \i(2).ip is given by 

K(2),^{n - rii) = \ph,^{n~m) + 5r,inli{2),v (36) 



with the electron-phonon coupling function \ph.ip{ri). 
Since EPI and E^_^ in Eq. (|34I35|) has a local form as 
a function of the angle <p) due to FSP, then the equation 
for Cjn,Lp has also a local form. This means that differ- 
ent points on the Fermi surface are decoupled. In that 
case Cjn,ip depends on the local value (on the Fermi sur- 
face) of the gap l^n,ip ~ Aq cos 2ip. This property alone 
is important in solving the ARPES non-shift puzzle. So, 
at the nodal point (tp = 7r/4) one has A„^i^ — and the 
quasiparticle spectrum given by 



E-^^-^k{E,Ar. 



0) = 



(37) 



is unaffected by superconductivity, i.e. the kink is not- 
shifted in the superconducting state, as shown in Fig.|H^. 
This is exactly what is seen in the experiment of the Shen 
group [23, - see Fig.|31 The FSP theory also predicts: 

(1) that E2(= ImS) has a knee-like shape for uj < uj''^^, 
see Fig.Eb, exactly as seen recently by ARPES |23|, |6a |: 

(2) that for uj > Wp^^"'' the EPI contribution to S2 is 
constant while its slope in this region is determined by 
the Coulomb scattering giving a small coupling constant 
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FIG. 6: Nodal direction: (a) The quasiparticle-spectrum 
"^(Ck) s-iid (b) the imaginary self-energy JmE(^ = 0,uj) in the 
nodal direction {ip = 7r/4) in the superconducting (T = 0.2 
meV) and normal {T — 6 meV) state, fic = 400 meV is the 
cutoff in Sc. After Ref. 



FIG. 7: Anti-nodal point: (a) The quasiparticle-spectrum 
ci;(5j^) and (b) the imaginary self-energy ImT,{^ = 0,^) in 
the anti-nodal direction {(p = 0; n/2) in the superconducting 
|7^ = 0.2 meV) and normal (T = 6 meV) state. After Ref. 



Ac < 0.4. Such a behavior is just what has been ob- 
served in ARPES spectra |23], |62ji thus confirming the 
correctness of the FSP theory. Let us point out that the 
experimental ARPES curves for Si and E2 are smoother 
than (our) theoretical curves, since we study the problem 
in the very simplified model with the slightly broadened 
Einstein spectrum. Since the tunnelling and ARPES ex- 
periments indicate that the broad spectrum of phonons 
contribute to EPI, then the quantitative analysis must 
rely on the realistic Eliashberg spectral function which is 
at present unknown. 

In the case of antinodal point {ip ~ 7r/2) the calcula- 
tions show, that there is a singularity at 40 meV (Egmg) 
in the quasiparticle spectrum in the normal state - see 
Fig. El The analytic and numeric calculations of Eq. H34|l 
show that this singularity is shifted by Aq in the super- 
conducting state, i.e. Esing ^ Esing + Ao- This Is cxactly 
what is seen in the recent experiment on BSCCO - see 
Fig. 121 where the singularity of the normal state spectrum 
at 40 meV is shifted to (65 — 70) meV in the supercon- 
ducting state, since Aq « (25 — 30) meV. We stress that 
the couphng constants Xph (> 1) and Ac (~ 0.4) were es- 
timated from Re S by knowing that Re S = —(Xph+Xc)LO 
for u; < w")^'^ and Re S = -Ac^^ for w")^'^ < w < (the 



band width) . The FSP theory explains in the natural way 
also the peak- dip-hump structure in ^(k , uj) if Xph > 1, 
as observed experimentally- see more in j37|. 

In conclusion, the different behavior of the nodal and 
the anti-nodal kinks in the superconducting state can be 
consistently described by the existence of the forward 
scattering peak in EPI and other charge scattering pro- 
cesses. We stress that the FSP theory is up to now the 
only one that succeeded in resolving the non-shift. 



2. ARPES isotope effect 

Recent ARPES spectra in the optimally doped Bi- 
2212 taken near the nodal and the anti-nodal points |2^ 
show after oxygen isotope substitution ^^O -^^^ O a pro- 
nounced isotope effect in the real part of the quasiparticle 
self-energy. This result, if confirmed, will be a crucial ev- 
idence for the importance of EPI in cupratcs. 

A. Isotope effect along the nodal direction. The 
ARPES results for the isotope effect in optimally doped 
Bi2212 samples are as follows [24] : {1) there is a kink 
in the quasiparticle spectrum at uJkjo — 70 meV which 
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is unshifted in the SC state; (2) there is a red shift 

SuJk,70 (10 - 15) meV of the kink for the ^^O O 

substitution; {3) the isotope shift of the self-energy is 
more pronounced at large energies lo — 100 — 300 meV. 

The results {l)-{2) can be explained by the EPI the- 
ory with FSP in a natural way. As it is already dis- 
cussed, the theory predicts that the kink in the super- 
conducting state is non-shifted with respect to that in 
the normal state, what is due to the FSP in theJEPI. Re- 
lated to the isotope effect, the theory predicts 67] that in 
the nodal direction after the substitution ^^O — O ( 
woie — ^kink ~ 70 meV) there is a red-shift of i?eE(k, cu) 
in the normal and superconducting state, due to the 
isotope dependence of the maximum of ReT,(k,Lu) at 
CO f» iOkink- The isotope effect is more pronounced at 
higher energies. This is shown in Fig. |S1 for the Debye 
model with the spectral function a'^F{uj) — X{uj/ujo)'^ 
with LOO — {k/MoY^"^ {Mo is the oxygen mass) - see the 
left panels in Fig. |S1 

The pronounced isotope effect is also present in the 
quasiparticle dispersion uj{k) with the read shift of the 
kink energy upon the substitution ^^O O. This is 
shown in the right panels of Fig. El The property (5) 
is difficult to explain by the Eliashbcrg-like theory we 
have used above. Some speculations that the pronounced 
isotope effect at high energies is of polaronic origin wait 
for further theoretical elaborations. 

B. Isotope effect near the anti-nodal point. The 
main experimental results for the ARPES isotope effect 
near the anti-nodal point k^w « (0, tt) are |2^: (i) there 
is a kink in the quasiparticle spectrum around u)k,io ~ 40 
meV which is shifted in the SC state by A(k^jv, w) « 30 
meV; (2) there is a red shift of the kink in the SC state 
with the energy change Suk ^ —5 meV at T = 25 K, 
which is smaller than the corresponding one in the nodal 
direction; (3) there is an inverse IE at higher energies 
io > UkAO, i.e. t^(efc;'^0) > c^o(^fc;'^0); U) in the 
normal state at T = 100 K there is practically no red 
shift for the ^^O O substitution. 

Theoretical calculations of the self-energy in the anti- 
nodal region have been done for the Debye spectrum '67] , 
only. Since the density of states near the anti-nodal 
point is higher than at the nodal point, we assume a 
larger (than in the nodal case) EPI coupling constant 
•^ep.kAjv = 2. In Fig. ^1 (top) it is seen that in the SC 
state (at T = 25 K) the kink in the quasiparticle energy 
is shifted by the gap A(k^Ar) and that there is a moderate 

red-shift ( (5-8) meV) of the kink for the ^^O -^^^ O 

substitution. Note that this value is only accidentally 
in agreement with the experimental one ( ~ — 5 meV). 
However, the theory predicts smaller (compared to exper- 
iments) inverse isotope effect {io{^f.;^^ O) > uJo{£,k]^^ O)) 
which occurs at higher energies. It seems that the ori- 
gin of this behavior lies in the smaller gap value A(k^Ar) 
of (dashed line in Fig. ^ than of ^^O (bold line). 
In spite of the fact that these results resemble qualita- 
tively the experimental results (i)-(5) there is quanti- 
tative difference between the theory and experiment at 
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FIG. 8: The shift of i?eE(k, ui) in the nodal direction upon 
the isotope substitution ^^O -^^^ O calculated using the De- 
bye model with loqib = 60 meV for various EPI coupling 
constants X at T ^ 25 K (top) and T = 100 7^ (bottom). 
From Ref. Igl. 



high energies, the origin of which is unknown at present. 

The FSP theory predicts that the isotope effect at the 
anti-nodal point is much less pronounced in the normal 
state (at 100 K) as shown in Fig. ^| These results 
resemble qualitatively the experimental results in ^2 j|. 
The origin of the small shift of the kink in the normal 
state lies in the lower antinodal kink energy and in the 
smearing effects at higher temperatures. 



3. Collapse of the anti-nodal elastic scattering in the 
superconducting state 

A number of ARPES experiments on optimally doped 
Bi-2212 give evidence for the significant elastic scattering 
which varies along the Fermi surface, being the smallest 
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FIG. 9: The quasiparticle energy uj{k — kp) in the nodal di- 
rection upon the isotope substitution ^^O O calculated 
using the Debye model with u>qi6 — 60 meV for various EPI 
coupling constants \ a,t T — 25 K (top) and T = 100 K 
(bottom). FromRef. 



along the nodal direction and the largest near the anti- 
nodal point, see |2lj | and references therein. The exper- 
iments indicate that the quasiparticle spectral function 
A(k, to) at the anti-nodal point is broad in the normal 
state and strongly sharpens in the superconducting state. 
On the other side, at the nodal point A(k, uj) has a stan- 
dard Lorentzian form which is almost unchanged in the 
superconducting state. 

This dramatic sharpening of the spectral function near 
the anti-nodal point (7r,0) at T < Tc can be explained 
by assuming that there is a forward scattering peak in 
the elastic impurity scattering as it was done in [s^- It 
is easy to see from Eq. (|34I35II . which describe the ex- 
treme forward scattering, that in the case of the FSP 
in impurity potential one has 7i = 72- As a result. 
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FIG. 10: The quasiparticle energy uj{k — kF) in the anti-nodal 
region under the isotope substitution ^^O -^^^ O in the Debye 
model with cjQie = 40 meV for the EPI coupling constants 
Aep.k^jv = 2 at T=25 K (top) and T= 100 K (bottom), kp 
and VF are the Fermi momentum and velocity, respectively. 



one finds that the impurity scattering rate rimp(k, w) = 
r„(k, cj) + ra(k, a;) = for I w |< Aq for any kind of 
pairing (s- p- d-wave etc.) since the normal (r„) and the 
anomalous (Fq) scattering rates compensate each other 
in the superconducting state - the collapse of the elastic 
scattering rate. This result is a consequence of the An- 
derson theorem [G^-fi^- In such a case d-wave pairing 
is unaffected by impurities - there is a negligible reduc- 
tion in Tc, as was first elaborated in 68] and 69], and 
further studied in '75]. The physics behind this result is 
rather simple. Forward scattering means that electrons 
scatter into a very small region in the k-space, so that 
at the most part of the Fermi surface there is no mix- 
ing of states with different signs of the order parameter 
A(k), and the detrimental effect of impurities is reduced. 
Only for states near the nodal points there is mixing but 
since A(k) is a small there is only small reduction in 
Tc 33]. The collapse of the elastic scattering rate in the 
ARPES spectra was recently considered in |73, where a 
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number of interesting results are reported. 

The idea of the forward scattering peak in the non- 
magnetic impurity scattering was appHed recently in [tJI 
in the study of the anomalous temperature dependence 
of the Hall angle in optimally doped HTSC. By taking 
into account (besides the inelastic scattering) also the 
small- angle elastic scattering on non-magnetic impuri- 
ties in [7l| it was explained the experimental finding in 
Y Ba2Cu3-j:Znx07-s that cotOniT) = axx/c^xy ~ T"^ 
{Qh is the Hall angle, H is the applied magnetic field) 
for IQQK <T < 300K. This is contrary to the expected 
linear behavior (for cotQniT)), since in that temperature 
interval one has p{T) ~ T. 

In conclusion, in order to explain ARPES results in 
cuprates it is necessary to take into account: (1) EPI 
interaction, since it dominates in the quasiparticle scat- 
tering in the frequency region responsible for pairing in 
cuprates; (2) effects of elastic nonmagnetic impurities 
with FSP; (3) the Coulomb interaction which dominates 
at higher energies uj > ujph- In this respect, the pres- 
ence of ARPES kink and the knee-like shape of the spec- 
tral width are the smoking-gun experiments that strongly 
constraint possible theories. 

B. Nonmagnetic impurities and robustness of 
d-wave pairing 

In the presence of strong correlations the impurity po- 
tential is also renormalized, as it is mentioned above, and 
the effective potential in the Born approximation is given 
by M?(q) = 7c(A:F,q)u.?,o(q)' where Ui,o(q) is the single 
impurity scattering potential in the absence of strong cor- 
relations. This was proved first in [s^ and elaborated 
quantitatively in [b^ and [b^. Since the charge vertex 
7c(PF;q) is peaked at q = the potential Wi(q) is also 
peaked at q = 0. This means that the scattering ampli- 
tude contains not only the s-channel (as usually assumed 
in studying impurity effects in cuprates), but also the d- 
channel, etc. Based on this property the FSP theory suc- 
ceeded in explaining some experimental fa cts, such as: (i) 
suppression of the residual resistivity jsfl, UJ. This 
effect was observed in optimally doped YBCO, where 
the resistivity p{T) at T = K has a rather small value 
< 10 pVlcm.; (ii) robustness of d-wave pairing |6l|, 
in the presence of impurities and other defects in Cu02 
planes, such as Zn, Ga, O defects, etc., see [t^ . 

Early theories of the effect of nonmagnetic impurities 
on Tc in cuprates always assumed that Wi(q) = const, 
i.e. that only the s-wave scattering channel is present. 
Such a theory predicts that Tc{pi^c) = at a much smaller 

is) 

residual resistivity pi ^ ~ 50 /ificm, while the experimen- 
tal range is 200 ^l]cm< < 1500 p^cm The lat- 
ter experimental fact means that d-wave pairing in HTSC 
is much more robust than what the standard theory pre- 
dicts, and it is one of the smoking gun experiments in 
testing the concept of FSP in the charge scattering po- 
tential. It is worth mentioning that in a number of pa- 



pers the pair-breaking effect of non-magnetic impurities 
in HTSC was analyzed in terms of the impurity concen- 
tration Hi, i.e. the dependence Tdrii). However, rii is 
not the parameter which governs this pair-breaking ef- 
fect. The more appropriate parameter for discussing the 
robustness of d-wave pairing is the impurity scattering 
amplitude r{9, 9'), which can be related to the measured 
residual resistivity and to Tc{pi). The robustness of 
d-wave pairing in HTSC can be revealed only by study- 
ing the experimental curve Tc(pj), as was first recognized 
experimentally in [t^I and theoretically in [s^l, [68f . 

The robustness of d-wave pairing in HTSC was ex- 
plained first in |Q, where the FSP theory (S^, [sj was 
applied to this problem. We shall not go into details - 
see [g^, - but just give here the general formula for 
the Tc{pi) dependence in anisotropic (including uncon- 
ventional) superconductors. We assume that in an un- 
conventional (anisotropic) superconductor the order pa- 
rameter has the form A(0) — /S.oY{6) and generally one 
has {Y{9)) ^ {{Y*{e)Y{e)) = 1). In the case of FSP, 
the momentum dependent impurity scattering amplitude 

is given by T{e, 0') = r,(0, 9') + T dYd[9)Yd{9') + and 

the critical temperature Tc by 

\n^ = ^{\)-^{\ + {l-l3)x) 

led Z Z 



-{Y{9)f[m{-)~^>{x+-)]. (38) 

Here, Tco is the bare critical temperature in the absence 
of impurities, x = Ts/iTrTc, f3 = Td/Ts and {Y{9)) means 
averaging over the Fermi surface. If one takes into ac- 
count that the inelastic scattering, which is character- 
ized by the coupling constant Xz, screens the impurity 
scattering, then the parameter x in Eq. H38|l should be 
replaced by x/{l + Xz)- Note that Eq. (|^ holds in- 
dependently on the scattering strength Fg, F^, i.e. it 
holds in the Born limit as well as in the unitary limit. 
The residual resistivity p^ can be related to the trans- 
port scattering rate by p^ — AnTtr/i^pi, while the s-wave 
amplitude is related to Ttr by F^ = pTtr- The parameter 
p > 1 can be obtained from the microscopic model (in the 
t-J model p w 2 — 3) or can be treated as a fitting parame- 
ter, see more in 0- In the case of an unconventional pair- 
ing in the tetragonal lattice one has {Y{9)) = and the 
last term drops. For the s-scattering only (F^ = 0) one 
has l3 — and Tc{pj) should be suppressed very strongly 
contrary to the experimental results W^, see Fig. 1111 

The FSP theory of impurity scattering in the t-J model 
[g^ gives that the s-channel and d-channel almost equally 
contribute to the impurity scattering amplitude, since 
13 « 0.75-0.85 for doping 5 « 0.1-0.2. The dependence 
of I3{5) has been calculated for the t-J model, see [68l| . 

Since the d-channel in scattering is not detrimental 
for d-wave pairing the FSP theory predicts that TdPi) 

vanishes at much larger p^^, i-e. pfj^^^ » p^il, in good 
agreement with experiments, as shown in Fig. 1111 In 
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FIG. 11: The critical temperature Tc [K] of a d-wave su- 
perconductor as a function of the experimental parame- 
ter p^/ac[K], where is the residual resistivity and a = 
Stt'^ Xtr/uJpi. The case /3 = corresponds to the pre- 
diction of the standard d-wave theory with isotropic scat- 
tering (7^ . The experimental data ff^ are indicated 
by crosses for YBa2{Cui-xZnx)307-s, and by circles for 
Yi-yPiyBaiCuaOr^s. From Ref. 

the extreme forward scattering limit, Tcq would be un- 
renormalized, i.e. the Anderson theorem holds also for 
unconventional superconductors. In this case impurity 
scattering does not mix superconducting states with dif- 
ferent sites of A(k). 

Here, we have studied the effects of impurities and 
other defects placed in the Cu02 planes. The corre- 
sponding experiments can be explained by FSP in the 
impurity potential, which is due to strong correlations. 
In the case when impurities are placed out of these planes 
their effect is weaker and also there is a forward scatter- 
ing peak due to weak screening along the c-axis. The 
proposed theory in j69i| and Ea.f (P5|l holds for both 
cases equally. 



C. Nonadiabatic corrections of Tc 

Cuprates are characterized not only by strong correla- 
tions but also by the relatively small Fermi energy Ep^ 
which is not much larger than the characteristic (max- 
imal) phonon frequency wJJ)^^, i.e. Ep — 0.1 — 0.3 eV , 
^max ^ gQ xhe situation is even more pronounced 

in fullerene compounds A^Cqq, with Tc = 20 — 35 K, 
where Ep ^ 0.2 eV and wjjjf' ~ 0.16 eV. This fact im- 
plies a possible breakdown of the Migdal's theorem [t^ , 
which asserts that the relevant vertex corrections due to 
the E — P interaction are small if (lod/Ef) ^ 1. In that 
respect a comparison of the intercalated graphite KC% 
and the fullerene A^Cqq compounds, given in ■76'|, is very 



instructive, because both compounds have a number of 
similar properties. However, the main difference in these 
systems lies in the ratio uj^jEp, since {ujd/Ep) <^ 1 in 
KCg, while it is rather large (uJo/Ep) ~ 1 in A3 Ceo ■ 
Since the ratio uJu/Ep is not negligible in the fullerene 
compounds and in cuprates it is necessary to correct 
the Migdal-Eliashberg theory by vertex corrections due 
to EPI. It is known that these vertex corrections lower 
Tc in systems with isotropic EPI. However, in systems 
with FSP and with the cut-off qc « kp in the pairing 
potential these vertex corrections cause an appreciable 
increase in Tc. The calculations by the Pietronero group 
[76| provided two important results: (1) there is a dras- 
tic increase of Tc by lowering Qc = Qc/^kp, for instance 
Tc(Oc = 0.1) « 4Tc(Qc = 1); (2) Even small values of 
A < 1 can give large Tc. The latter results open a new 
possibility in reaching high Tc in systems with apprecia- 
ble ratio LOo/Ep and with FSP. The difference between 
the Migdal-Eliashberg and non-Migdal theories can be 
explained qualitatively in the framework of an approx- 
imative McMillan formula for Tc (for not too large A) 
which reads Tc « (tj)e"i/[^"^*l . The Migdal-EUashberg 
theory predicts Xme ~ A/(l + A) while the non-Migdal 
theory Q| gives XnM ~ A(l -I- A). For instance Tc ~ 100 
K in cuprates can be explained by the Migdal-Eliashberg 
theory if A ^ 2, while in the non-Migdal theory much 
smaller coupling constant is needed, i.e. A ~ 0.5. 

VI. PSEUDOGAP AND THE 
ELECTRON-PHONON INTERACTION 

The pseudogap (PG) problem is a very intriguing one 
and it is not surprising that a variety of theoretical mod- 
els have been proposed to explain it. We are not going 
to discuss all these in detail but only mention some that 
may have experimental support. The first one is based 
on the assumption that the PG phase represents the state 
with pre-formed pairs |77|. where the true critical tem- 
perature Tc is smaller than the mean- field one T*^^. In 
the region Tc < T < pre-formed pairs give rise 

to a dip in the density of states N{uj). This approach 
is physically plausible having in mind that cuprates are 
characterized by a short coherence length, quasi-two di- 
mensionality and proximity to the Mott-Hubbard insu- 
lating state. As we already discussed the latter gives 
rise to small phase stiffness K'^ of the superconduct- 
ing order parameter A =| A | e*'''. In 2D systems the 
Berezinskii-Kosterlitz-Thouless transition temperature, 
Tc = {'^/'2)Ks{Tc) S for small doping), can be much 
smaller than Tc ^^ ■ From the experimental side there is 
some support for this approach, at least in not very broad 
temperature region close to Tc. For instance, the spe- 
cific heat measurements [t^ (with smaller jumps at Tc) 
point to a non-mean field character of the superconduct- 
ing phase transition, particularly for underdoped com- 
pounds. As we have already mentioned in the Introduc- 
tion, ARPES measurements show that for Tc < T < T* 
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PG is d-wave like Apg(k) « Apg.o(cos fc^;— cos ky), like the 
superconducting gap, and that A^g q increases at lower 
doping. 

The second approach assumes that PG is due to exis- 
tence of an additional competing order, but usually with- 
out the true long-range order. For instance, the "spin- 
density wave" alias for strong antiferromagnetic fluctu- 
ations can also produce dip in N{uj) 79]. Related ap- 
proaches are based on RVB, orbital currents, d-wave or- 
der, etc., which are not discussed here. 

We would like to pay attention to two possible effects, 
which are due to FSP, that can also give rise to PG in the 
density of states. The first is related to EPI interaction 
with FSP that can produce a dip in N{llj) for ut < ft, 
where fi is the characteristic phonon frequency. The 
second is related to a novel type of fluctuations - internal 
Cooper pair fluctuations, due to the long-range pairing 
forces (alias for FSP) - due to EPI with FSP. 



2.5 
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1. Pseudogap due to phonons 

The EPI, in the Einstein model with the phonon fre- 
quency ri, was also studied in the extreme limit of FSP 
|8Q], where in leading order the spectral function is sin- 
gular, i.e. a^i^(k, k', (jj) ~ (5(k — 'k')5{Lo — f2). In such a 
singular case besides renormalization of frequency {lo —> 
Z(\i,Lu)uj) there is also significant renormalization of the 
quasiparticle energy (^j. xO^^'^)^k) the Eliashberg 
equations, giving a v ery interesting behavior of the den- 
sity of states N{uj) |80j- First, there is an increase of 

N{U) = 0)(> Nbarei^ = 0)) at W = 0, where Nbare{^^ = 0) 

is the bare density of states in the absence of EPI, see 
Fig. 112b . This fact is mainly due to renormalization of 
the quasiparticle energy (^j. — > xO^^'^)^k)- 

In such a case there is a pseudogap feature in the re- 
gion < w < f2, where N{uj) < Nbarei^^)- This 
pseudogap feature disappears at temperature compara- 
ble to the phonon energy fl. Note that the usual isotropic 
and short-ranged EPI does not renormalize the density 
of states in the normal state, i.e. N{uj) — Nbare{^)- 

It is interesting that without renormalization of the 
quasiparticle energy (^j^ ^j., i.e. x(k,w) — 1) there 
is also pseudogap behavior. The case x = 1 is probably 
realized in systems with finite width of FSP {5k <ti kp) 
and with particle-hole symmetry. In such a case one has 
N{uj = 0) < Nbare{'-^ — 0) and the pseudogap exist for 
LO <Vl, as it is seen in Fig. 112b. 

Let us mention that for x(k,w) ^ 1 the transport prop- 
erties are very peculiar due to the pseudogap behavior of 
N{lo) and the existence of peak at N{u) = 0). For in- 
stance, the resistivity p{T) is linear in T starting at very 
low temperatures, i.e. p{T) ~ T for (fi/SO) < T and this 
linearity extends up to several f2. The dynamic conduc- 
tivity i7i(w) shows the (extended) Drude-like behavior 
with the Drude width Ttr ~ T , for lo < T. The above 
numbered properties are in a qualitative agreement with 
experimental results in cuprates - see more in |8C| . 0- 
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FIG. 12: (a) The density of states iV(aj) in the FSP model for 
EPI with the dimensionless coupling l{= Vep/ttO.) = 0.1 for 
various t{= nT/Q) - case ~ x(k,i^)Ck- (b) N(u) for the 
same parameter 1 but for S^y^^ ~ ^k- From Ref. [80||. 



In this extreme FSP limit one can also calculate the 
mean- field critical temperature T^q^- In leading order 
of n/T^^^ (» 1) one has T^^^ « N{0)Vep = AiV(0)/4, 
where A = N{Q)Vep- In this case the maximal super- 
conducting gap is given by Aq = 2T^^^^ which is reached 
on the Fermi surface, while away from it the gap de- 
creases, i.e. A/j = Ao-\/l — The expression for 
T^Q^ tells us that it can be large even for A < 0.1, since 
in cuprates the bare density of states is Nf,areiO) ^ 1 
states /eV. It is apparent that in this order there is 
no isotope effect in Tjjf^ since a = 0. However, such 
an extreme limit is never realized in nature, but for a 
qualitative understanding of the self-energy it is a good 
starting point, since the effects of the finite width (fcc) 
of Q;^i^(k, k', w), whenever kn<€i kp, change mostly the 
quantitative picture - see |8(]|. In case of finite width 
of FSP when kcVp ^ the reduction of Tc is given by 
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rj^MF ^ T^F{i _ 7C(3)fcet;f /47r2T^^^). There is also re- 
duction of T*^^ due to retardation effects giving rise to fi- 
nite isotope effect a ^ 0. Such a singular model has 
repercussions on other properties. For, instance interest- 
ing calculations within more realistic FSP theory with 
the finite width fcc, but kc <C fcF, were done in ■81], where 
the FSP theory for EPl and the SFl theory (based on 
spin-fluctuation mechanism of pairing) were compared. 
It was shown in this interesting paper that the FSP the- 
ory can explain the appreciable increase of anisotropy 
gap-ratio R = A(7r, 0)/A(7r/2, 7r/2) when T T,., while 
the SFI theory fails. Furthermore, the FSP theory of 
EPI can explain the pronounced effect of orthorhombicity 
(a 7^ b) in YBCO on the gap ratio Aa/A;,, the anisotropy 
of penetration depth ^a/^l the supercurrent ratio in 
the c-axis Pb-YBCO junction. On the other hand, the 
SFI theory is ineffective here, since it predicts at least 
one order of magnitude smaller effects [sj , • 

2. Pseudogap due to internal Cooper pair fluctuations 



long-range pairing potential there are additional fluctu- 
ations of the internal degrees of freedom (r) . In the fol- 
lowing we sketch the analysis given in [83|. When the 
range of pairing potential is large, i.e. Tc > ^ (the super- 
conducting coherence length), fluctuations of the internal 
Cooper pair wave-function are important since they give 
rise to a large reduction of the mean- field quantities. In 
order to understand the physics of internal wave-function 
fluctuations we study a much simpler Hamiltonian - the 
so called reduced BCS Hamiltonian, 

HbCS = X^^kcLck^ - ^k-k'4T"^-kiC-k'iCk'i- 
ko- k,k' 

(40) 

Since we shall study excitations around the ground state 
we assume that there are no unpaired electrons which al- 
lows us to study the problem in the pseudo-spin Hamil- 
tonian (see |82| and references therein) 

Hps = Yl - ^ E ^k-k' {S+Si^, + S+Si^) 

ko- k,k' 



As it is stated above, the FSP theory also predicts 
existence of a long-range pairing force due to the renor- 
malization of EPI by strong correlations. This opens an 
additional possibility for PG. Due to FSP in EPI, the 
mean-field critical temperature has a non-BCS depen- 
dence, i.e. T^^^ = Vep/4: and it is inevitably reduced by 
phase and internal Cooper pair fluctuations. The latter 
are always present in systems with long-range attractive 
forces, i.e. with FSP, as it was pointed out first in 82]. It 
was shown there that such a long-ranged superconductor 
exhibits an additional class of fluctuations in which the 
internal structure of Cooper pair is soft. This leads to PG 
behavior in which the actual transition temperature Tc 
is greatly depressed from its mean- field value T^^^ . We 
stress that these fluctuations are not the standard phase 
fluctuations in superconductors. Because of the impor- 
tance, as well as due to internal beauty, of the approach 
given in "82"] we summarize briefly its main results. 

In the following the weak coupling limit for the pairing 
Hamiltonian is assumed 

dxV.I(x)e(p)^,(x) 



- j dxdxV(x-x')V'I(x)?/'l(x')?/'|(x')V'TW. (39) 

In the mean-field approximation (MFA) the or- 
der parameter A(x, x') is given by A(x, x') = 
F(x — x')('0j^(x')'0|(x)), which depends in fact on the 
internal coordinate r = x — x' and the center of mass 
R = (x + x')/2, i.e. A(x,x') = A(r,R). In standard su- 
perconductors (LTSC) with a short-range pairing poten- 
tial Vsr{^ - x') K. Vb(5(x - x') one has A(r,R) = A(R) 
and in that case only the spatial (R) fluctuations of the 
order parameter might be pronounced. In the case of 



= E 2ek5L - E ^k-k' {SlSl, + SI, SI), (41) 

ko- k,k' 

where the pseudo-spin 1/2 operators S^^, S^^ = {S^^Y 
are given by 

= ^(c^TCkT - clkiC-kj - 1), 

^k. = 4TSiki- (42) 

We see that Eq. (|41|l belongs to the class of Heisen- 
berg ferromagnetic Hamiltonians (since Vk_k' > 0) with 
the lattice in the Brillouin zone. It is well-known for 
the Heisenberg model that in the case of long-range 
forces Vk-ki = const in a large part of the k-space, then 
the phase fluctuations are negligible and thermodynamic 
properties are well described in MFA 

i?A/FA = -"E^kSk (43) 
k 

with the mean-field hk given by 

hk - -2ekZ+ E ^k-k' (^^'X+^^,y) , (44) 

k' 

where x, y and z are unit vectors. Since x- and y- 
axis are equivalent one can look for hk in the form 
hk = — 2(^i(.z-t-2AkX, where the order parameter Ak is 
the solution of equation 

Ak = E^k-k^(^k')-E%^tanh^, (45) 



with i?k — Y^k + ^k- stress that the long-range 
force (Vk-k' — const) in the k-space produces in the 
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real space short-range BCS-like force Vbcs{'^ — x') « 
Vo<5(x — x'). This situation is realized in LTSC supercon- 
ductors where the metallic screening makes the pairing 
potential short-ranged. 

For the long-range attractive forces the function Vk_k' 
is sharply peaked at | k — k' |= 0, for instance in the 
extreme FSP case (see Section 6.2) one has Vk-k' = 
Vb(5(k — k'). In such a case we deal with the short-range 
Heisenberg model in the k-space. It is well known that in 
such a system low-laying excitations, such as the Gold- 
stone magnons, tend to destroy the long-range order, i.e. 
in these systems one has Tc ^ T*^^. In fact these fluctu- 
ations correspond to the internal Cooper pair fluctuations 

In the following we analyze s-wave pairing only, where 
the solution of Eq. ^ gives T*^^ = Vb/4 and Ao = 
2T^'^. (Note, that in the BCS case one has Aq = 1.76 
T^^^ .). The coherence length is defined by ^ = vp/ttAq. 
As we said, the important fact is that in the case of long- 
ranged superconductors the Heisenbcrg-likc Hamiltonian 
in the momentum space is short-ranged giving rise to low- 
laying spin-wave spectrum. The latter spectrum are in 
fact the low-energy bound states (excitons) which loosely 
correspond to the low-energy collective modes (in the 
true many-body theory based on Eq. (|^ ). This prob- 
lem was studied in [221 for the long-range (but finite) 
potential V{r) — Vb exp{— r^/2r^} (its Fourier trans- 
form is Vfc — (27rr^)Vb exp{— fc^r^/2}), and it was found 



Bu = 



that a large number Nc> 



Trkprc/G^ (for rc ^ of 



the excitonic like collective modes wj^j^ exist at zero mo- 
mentum. These excitonic modes lie between the ground 
state and the two-particle continuum, i.e. their energies 
are a; < 2Ao. Note that since we assume that Aq ^ Ep, 
the system is far from the Bose-Einstein condensation 
limit. 

The above analysis in terms of the pseudo-spins is use- 
ful for physical understanding, but the full many-body 
fiuctuation problem, which is based on the Hamiltonian 
in Eq. (|39|l . is studied in jsj where the Ginzburg- 
Landau (G-L) equation is derived for the long-ranged 
superconductor. Due to fiuctuations of the internal wave- 
function the G-L free-energy functional F{A(R, k)} for 
the order parameter A(R, k) = J (irA(R — r/2,R + 
r/2) exp{— ikr} has a more complicated form 

F{A(R, k)} = ^ /" dR{Ak I A(R, k) p \ A(R, k) 



^|akA(R,k) P+^|aKA(R,k) PI, 



(46) 



where M = r'^Vo and 



Vn 



tanh(/3gj2) 



1 _ /32«|sinh(/34/2) 



tanh(/3Ck/2) 
S^k 



/3 



16C^cosh^(/3^k/2) 



(47) 



The term with the partial derivative 9k is a direct 
consequence of the long-ranged pairing potential, and it 
describes fluctuations of the internal Cooper pair wave- 
function. The effect of these fluctuations, described by 
the free-energy functional in Eq. H46|) , was studied in the 
Hartree-Fock approximation in the limit Tc ^ ^, and a 
large reduction of the mean-field critical temperature is 
found 



MF 



irc/0' 



(48) 



2TOk 32^1, cosh'' (/3ek/2) 



The latter result means that Tc in the long-range su- 
perconductors is controlled by thermal fluctuations of col- 
lective modes which is in contrast with the 3D short-range 
(BCS-like) superconductivity, where phase fluctuations 
dominate but only slightly suppress T*^^. In the tem- 
perature interval Tc < T < Tj^^^, the system with long- 
range pairing force is in PG regime where the electrons 
are paired but there is no long-range phase coherence. 
The latter sets in only at T < Tc. We shall not further 
discuss this interesting approach but only stress that it 
can be generalized by including repulsive interactions, for 
instance due to spin fluctuations and residual Coulomb 
interaction. This will be studied elsewhere [83l |. 

In conclusion, FSP in EPI gives rise to the long-range 
pairing (in real space) which produces internal fluctua- 
tions of Cooper pairs. These are soft excitonic modes 
of the internal Cooper wave function which may reduce 
Tc strongly. In the region Tc < T < T^^ PG phase 
is realized. In this approach PG possess the same sym- 
metry as the superconducting gap. Since the cuprates 
are systems near the Mott-Hubbard insulating state due 
to strong correlations, this is responsible for very small 
phase stiffness and strong phase fluctuations. Addition- 
ally, strong correlations (in conjunction with the long- 
range Madelung interaction due to ionic-metallic struc- 
ture of cuprates) make EPI long-ranged, which causes 
strong internal fluctuations of Cooper pairs. In a realis- 
tic theory both type of fluctuations must be taken into 
account 1831. 



VII. SUMMARY AND CONCLUSIONS 

A number of experiments, such as optics {IR and Ra- 
man), transport, tunnelling, ARPES, neutron scatter- 
ing, etc. give convincing evidence that EPI in cuprates 
is sufficiently strong and contributes to pairing. The 
ARPES experiments made a breakthrough in the physics 
of cuprates, since they allowed a direct study of quasi- 
particle properties, such as self-energy effects. These ex- 
periments also give evidence for the presence of strong 
correlations which modify EPI not only quantitatively 
but also qualitatively. The most spectacular result of 
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the EPI theory in strongly correlated systems is the ap- 
pearance of FSP in EPI, as well as in other charge scat- 
tering processes such as the residual Coulomb interac- 
tion, scattering on non-magnetic impurities - the FSP 
theory [H, 0. FSP is especially pronounced 

at lower doping 6. This fundamental result allows us 
to resolve at least qualitatively a number of experimen- 
tal facts which can not be explained by the old the- 
ory based on isotropic Migdal-Eliashberg equations for 
EPI. The most important predictions of the FSP theory 
of EPI and other charge scattering processes are: (1) 
the transport coupling constant Xtr (entering the resis- 
tivity, g ^ XtrT) is much smaller than the pairing A, i.e. 
Atr < A/3; (2) the strength of pairing in cuprates is ba- 
sically due to EPI, while the residual Coulomb repulsion 
(including spin-fluctuations) triggers d-wave pairing; (3) 
d-wave pairing is very robust in presence of non-magnetic 
impurities with FSP; (4) the nodal kink in the quasiparti- 
cle spectrum is not-shifted in the superconducting state, 
while the anti-nodal singularity is shifted by the maxi- 
mal superconducting gap; (5) there is a collapse of the 
elastic impurity scattering rate in the superconducting 
state at the anti-nodal point giving rise to sharpening of 
the ARPES features; (6) internal fluctuations of Cooper 
pairs, due to the long-ranged scattering potential, can 
additionally increase the temperature range of the PG 
effect. 

We point out the following two facts that stem from the 
theoretical analysis: (i) FSP in EPI of strongly correlated 
systems is a general phenomenon and it affects electron 
coupling to all phonons; (ii) the existence of FSP in EPI 
is confirmed numerically by Monte Carlo calculations for 
the Hubbard- Holstein model with finite U by exact 
diagonalization j^^l, as well as by some other methods 
|6^. Thanks to these numerical results which confirm 
the correctness of the FSP theory [30l|. [3ll|. Jgl. [rf. a 
number of new results were recently reported '84'| , which 
are in favor of strong EPI interaction and the FSP theory. 



We would like to stress that the calculation of the mo- 
mentum dependent EPI coupling constant Ak in cuprates 
is a big challenge for the theory. The latter must in- 
clude properly: (a) ionic-metallic coupling due to the 
long-range Madelung energy as well as covalent coupling, 
and (b) strong electronic correlations. From this point of 
view, one can say that we still do not have a quantitative 
microscopic theory of pairing in cuprates. 

In the last several years a large number of papers has 
been devoted to the study of spin-fluctuation interaction 
as a mechanism of pairing in cuprates. In spite of much 
effort led by the greatest authorities in the field, which 
have opened some new research directions in the theory 
of electron magnetism, there is no theoretical evidence for 
the effectiveness of non-phononic mechanism of pairing. 
Moreover, recent numerical calculations 13 show that 
probably there is no high-temperature superconductivity 
in the t- J model. If it exists at all, its Tc is extremely low. 
Finally, numerical calculations in the Hubbard model 
show that the repulsive Hubbard interaction is unfavor- 
able for high-temperature superconductivity, contrary to 
the attractive interaction which favors it. 

In conclusion, we would like to stress again: in order 
to explain the high critical temperature in cuprates, it is 
necessary to include EPI which is renormalized by strong 
electronic correlations. We hope that in the future such 
projects will be supported much more. 
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